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Abstract 

We propose the Hyper-Kamiokande (Hyper-K) detector as a next generation underground water 
Cherenkov detector. It will serve as a far detector of a long baseline neutrino oscillation experiment en- 
visioned for the upgraded J-PARC, and as a detector capable of observing - far beyond the sensitivity 
of the Super-Kaniiokande (Super-K) detector - proton decays, atmospheric neutrinos, and neutrinos from 
astronomical origins. The baseline design of Hyper-K is based on the highly successful Super-K, taking full 
advantage of a well-proven technology. 

Hyper-K consists of two cylindrical tanks lying side-by-side, the outer dimensions of each tank being 
48 (W) X 54 (H) X 250 (L) m^. The total (fiducial) mass of the detector is 0.99 (0.56) miUion metric tons, 
which is about 20 (25) times larger than that of Super-K. A proposed location for Hyper-K is about 8 
km south of Super-K (and 295 km away from J-PARC) at an underground depth of 1,750 meters water 
equivalent (m.w.e.). The inner detector region of the Hyper-K detector is viewed by 99,000 20- inch PMTs, 
corresponding to the PMT density of 20% photo-cathode coverage (one half of that of Super-K) . 

Hyper-K presents unprecedented potential for precision measurements of neutrino oscillation parameters 
and discovery reach for CP violation in the lepton sector. With a total exposure of 5 years (one year 
being equal to 10^ sec) to a 2.5-degree off-axis neutrino beam produced by the 1.66 MW J-PARC proton 
synchrotron, it is expected that the CP phase 5 can be determined to better than 18 degrees for all possible 
values of 5 and CP violation can be established with a statistical significance of 3cr for 74% of the 5 parameter 
space if sin^ 20i3 > 0.03 and the mass hierarchy is known. If sin'^ 20i3 is as large as 0.1 the mass hierarchy 
can be determined with more than 3cr statistical significance for 46% of the 5 parameter space. In addition, 
a high statistics data sample of atmospheric neutrinos will allow us to extract the information on the mass 
hierarchy and the octant of 6*23. With a full 10 year duration of data taking, the significance for the mass 
hierarchy determination is expected to reach Sct or greater if sin^ 62^ > 0.4. 

Hyper-K can extend the sensitivity to nucleon decays beyond what was achieved by Super-K by an order 
of magnitude or more. The sensitivities to the partial lifetime of protons for the decay modes of p — > e+Tr*^ 
and p — > PK^ are expected to exceed 1 x 10"^^ years and 2 x 10'^'' years, respectively. This is the only known, 
realistic detector option capable of reaching such a sensitivity for the p e^n'^ mode. 

The scope of studies at Hyper-K also covers high precision measurements of solar neutrinos, observation of 
both supernova burst neutrinos and supernova relic neutrinos, dark matter searches, and possible detection 
of solar fiare neutrinos. The prospects for neutrino geophysics using Hyper-K are also mentioned. 



|mcisato@suketto. icrr.u-tokyo.ac.jp 
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FIG. 1. Schematic view of the Hyper-Kamiokande detector. 
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The Hyper-Kamiokande (Hyper-K or HK) detector is proposed as a next generation under- 
ground water Cherenkov detector that serves as a far detector of a long baseUne neutrino oscillation 
experiment for the J-PARC neutrino beam and as a detector capable of observing proton decays, 
atmospheric and solar neutrinos, and neutrinos from other astrophysical origins. The baseline de- 
sign of Hyper-K is based on the well-proven technologies employed and tested at Super-Kamiokande 
(Super-K or SK). Hyper-K consists of two cylindrical tanks lying side-by-side, the outer dimensions 
of each tank being 48 (W) x 54 (H) x 250 (L) m'^. The total (fiducial) mass of the detector is 0.99 
(0.56) million metric tons, which is about 20 (25) times larger than that of Super-K. A proposed 
location for Hyper-K is about 8 km south of Super-K (and 295 km away from J-PARC) and 1,750 
meters water equivalent (or 648 m of rock) deep. The inner detector region is viewed by 99,000 
20-inch PMTs, corresponding to the PMT density of 20% photo-cathode coverage (one half of 
that of Super-K). The schematic view of the Hyper-K detector is illustrated in Fig. [ij Table [l] 
summarizes the baseline design parameters of the Hyper-K detector. 



TABLE I. Detector parameters of the baseline design. 



Detector type 




Ring-imaging water Cherenkov detector 


Candidate site 


Address 


Tochibora mine 






Kamioka town, Gifu, JAPAN 




Lat. 


36°21'08.928"N 




Long. 


137°18'49.688"E 




Alt. 


508 m 




Overburden 


648 m rock (1,750 m water equivalent) 




Cosmic Ray Muon flux 


1.0 - 2.3 X 10"^ sec-icm-2 




Off-axis angle for the J-PARC v 


2.5° (same as Super-Kamiokande) 




Distance from the J-PARC 


295 km (same as Super-Kamiokande) 


Detector geometry 


Total Volume 


0.99 Megaton 




Inner Volume (Fiducial Volume) 


0.74 (0.56) Megaton 




Outer Volume 


0.2 Megaton 


Photo-multiplier Tubes Inner detector 


99,000 20-inch PMTs 






20% photo-coverage 




Outer detector 


25,000 8-inch </> PMTs 


Water quality 


light attenuation length 


> 100 m @ 400 nm 




Rn concentration 


< 1 mBq/m'^ 
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Hyper-K provides rich neutrino physics programs. In particular, it has unprecedented potential 
for precision measurements of neutrino oscillation parameters and discovery reach for CP violation 
in the lepton sector. With a total exposure of 5 years (1 year being equal to 10'' sec) to a 2.5-degree 
off-axis neutrino beam produced by the 1.66 MW J-PARC proton synchrotron, it is expected that 
the CP phase 6 can be determined to better than 18 degrees for all values of 6 and that CP 
violation can be established with a statistical significance of 3a for 74% of the 5 parameter space if 
sin^ 2^13 > 0.03 and the mass hierarchy is known. If the mass hierarchy is unknown, the sensitivity 
to the CP violation is somewhat reduced due to degeneracy. For a large value of sin^ 2^13, it is also 
possible to determine the mass hierarchy for some of 6 with this program alone. If sin^ 2^13 = 0.1, 
the mass hierarchy can be determined with more than 3a statistical significance for 46% of the S 
parameter space. The recent result of sin^ 2^13 > 0.03 obtained by the T2K experiment [1] boosts 
the expectation of discovery of CP violation by Hyper-K. 

The high statistics data sample of atmospheric neutrinos obtained by Hyper-K will allow us to 
extract information on the mass hierarchy and the octant of ^23- With a full 10 year period of data 
taking, the significance for the mass hierarchy determination is expected to reach 3a or greater if 
sin^ 2^13 > 0.04 and sin^ ^23 > 0.4. If sin^ 2^23 is less than 0.99, it is possible to identify the octant 
of 023, i.e. discriminate sin^ ^23 < 0.5 from > 0.5. 

Hyper-K extends the sensitivity to nucleon decays far beyond that of Super-K. The sensitivity 
to the partial lifetime of protons for the decay mode p — )• e~^7r^, the mode considered to be most 
model independent, is expected to be 1.3 x 10^^ years. It is the only realistic detector option 
known today able to reach this sensitivity. The sensitivity for the decay mode p — t- i>K~^ , the mode 
favored by super symmetry (SUSY) models, reaches 2.5 x 10'^^ years, and therefore Hyper-K would 
discover proton decay if some of the SUSY models are correct. 

Hyper-Kamiokande functions as an astrophysical neutrino observatory. If a core collapse super- 
nova explosion occurs halfway across our galaxy, the Hyper-K detector would detect approximately 
170,000~260,000 neutrinos as a ~ 10 second long burst. This very large statistical sample should 
at last reveal the detailed mechanism of supernova explosions. For instance, the onset time of the 
explosion can be determined with an accuracy of 0.03 milliseconds, which is a key information to 
study the first physical process of the explosion (p + e~ — )• n + Ug), allowing examination of the 
infall of the core and the ability to see the precise moment when a new neutron star or black hole 
is born. The sharp risetime of the burst in Hyper-K can also be used to make a measurement 
of the absolute mass of neutrinos. Because of non-zero masses, their arrival times will depend on 
their energy. The resulting measurement of the absolute neutrino mass would have a sensitivity 
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TABLE II. Physics targets and expected sensitivities of the Hyper-Kamiokande experiment, aso is the 
WIMP-proton spin dependent cross section. 



Pln sics Tar<>ot 



Sciisiti\i( \- 



Couditious 



Neutrino study w/ J-PARC v 

— CP phase precision < 18° 



— CPV 3c7 discovery coverage 



74% (55%) 
74% (63%) 



1.66 MW X 5 years (1 year = 10^ sec) 
@ s'^26is{= sin^ 26I13) > 0.03 and 
mass hierarchy (MH) is known 
@ 3^2^13 = 0.1, MH known(unknown) 
@ 5^2^13 = 0.03, MH known(unknown) 
(I; .s-26'i:i = 0.01. Mil kii()\vii(uukiunvu) 



Atmospheric neutrino study 

— MH determination > 3cr CL 

— ^23 octant determination > 90% CL 



10 years observation 

@ 0.4 < 5^6123 and 0.04 < s'^2exs 

@ 3^2923 < 0.99 and 0.04 < 8^2613. 



Nucleon Decay Searches 
- p e+ + 7r° 



10 years data 



1.3 X 10^5 yrs (90% CL) 
5.7 X 10^4 yrs (3a CL) 
2.5 X 10^4 yrs (90% CL) 
1.0 X 10^4 yrs (3a CL) 



Solar neutrinos 

-^Bv from Sun 

— *B 1/ day /night accuracy 



200 iv's / day 
< 1% 



7.0 MeV threshold (total energy) w/ osc. 
5 years, only stat. w/ SK-I BG x20 



Astrophysical objects 

— Supernova burst v 

— Supernova relic u 

— WIMP annihilation at Sun 



170,000~260,000 i/'s 

30~50 z/'s 

830 iv's / 10 years 

asD = 10-39cm2 
asD = lO-^Ocm^ 



@ Galactic center (10 kpc) 
@ M31 (Andromeda galaxy) 

5 years observation 

@ MwiMP = 10 GeV, XX bb dominant 

@ MwiMP = 100 GeV, XX ^ W+W~ dominant 



of 0.5 — 1.3 eV/c^, and does not depend on whether the neutrino is a Dirac or Majorana particle. 
Hyper-K is also capable of detecting supernova explosion neutrinos from galaxies outside of our 
own Milky Way; about 7,000-10,000 neutrinos from the Large Magellanic Cloud and 30-50 even 
from the Andromeda galaxy. 

Detection of supernova relic neutrinos (SRN) is of great interest because the history of heavy 
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element synthesis in the universe is encoded in the SRN energy spectrum. With gadohnium added 
to water, a neutron produced by the inverse beta process {ug + P ^ + n) , which is the predom- 
inant interaction mode for the SRNs, can be tagged by detecting gammas from the Gd(n,7s)Gd 
reaction. Doing so greatly reduces backgrounds and opens up the SRN energy window, improving 
the detector's response to this important signal. Our study shows that Hyper-K with 0.1% by mass 
of gadolinium dissolved in the water is able to detect as many as 830 SRNs in the energy range of 
10-30 MeV for 10 years of livetime. This large sample will enable us to explore the evolution of 
the universe. 

Dark matter can be searched for in Hyper-K as was done in Super-K. Neutrinos emitted by 
weakly interacting massive particles (WIMPs) annihilating in the Sun, Earth, and galactic halo 
can be detected using the upward-going muons observed in Hyper-K. A sensitivity to the WIMP- 
proton spin dependent cross section would reach 10~^^ (lO"^*') cm'^ for a WIMP mass of 10 (100) 
GeV and 5 years of livetime. 

Table [III summarizes the physics potential of Hyper-K. This document serves to define the scope 
of the Hyper-K project, describe a baseline design of the detector, and make a physics case for its 
construction. The required R&D items, cost, and schedule for constructing the Hyper-K detector 
will be provided in a separate documentation later. 

The Science Council of Japan announced the "Japanese Master Plan of Large Research Projects" 
for the first time on 17 March 2010 [2j. The plan includes 43 projects selected from 7 fields of 
science: Humanities and Social Sciences; Life Sciences; Energy, Environmental and Earth Sci- 
ences; Material and Analytical Sciences; Physical Sciences and Engineering; Space Sciences; and 
Information Sciences. There are a total of five large research facility projects listed in physical 
sciences. Among them is "Nucleon Decay and Neutrino Oscillation Experiments with Large Ad- 
vanced Detectors." This is a project to advance neutrino physics/astronomy and to search for 
nucleon decays using a large water Cherenkov detector that is approximately 20 times larger in 
volume than Super-Kamiokande and/or a large liquid argon detector. The Hyper-K detector is 
one of the leading options for this project. 

We are also well aware of other detector technology options such as liquid argon and of activities 
in other regions such as LBNE [3j in the US, and LAGUNA [4J in Europe. It is our firm intention 
to make the Hyper-K project completely open to the international community and contribute to 
the world-wide effort to make a strong neutrino physics program. 
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I. INTRODUCTION 
A. Physics case 

The goal of particle physics is to discover and understand the fundamental laws of nature. The 
Standard Model (SM), which is the current paradigm of elementary particles and their interactions, 
gives a successful account of the experimental data to date [5]. Yet, deeper insights are still needed 
to answer more profound questions. For instance, why does there exist a gauge structure of 
SU{3)c (d) SU{2)i (g) U{1)y among the strong, weak and electromagnetic interactions? Why is 
there a three generation structure of fundamental fermions and what are the origins of the masses 
and generation mixings of quarks and leptons? To address these questions physics beyond the SM 
(BSM) is required. 

The discovery of neutrino oscillations by the Super-Kamiokande (Super-K or SK) experiment in 
1998 [6] opened a new window to explore BSM physics. Evidence of neutrino oscillations is the only 
experimental proof known today that shows the existence of BSM physics at work. The mixing 
parameters of neutrinos, though not yet fully determined, were found to be remarkably different 
from those of quarks, which suggests the presence of an unknown flavor symmetry waiting to be 
explored. Extremely small masses of neutrinos compared with those of their charged partners lead 
to the preferred scenario of a seesaw mechanism, in which small neutrino masses are a reflection 
of the ultra-high energy scale of BSM physics. Furthermore, recent theoretical works point to the 
intriguing possibility that CP asymmetry originating from flavor mixing among the three genera- 
tions of neutrinos might have played an important role in creating the observed matter-antimatter 
asymmetry in the universe. Therefore, to explore the full picture of neutrino masses and mixings 
and to observe CP asymmetry in the neutrino sector are among the most important and urgent 
subjects in today's elementary particle physics world. The indication, obtained recently by the 
T2K experiment [1], that the mixing angle parameter between the first and third generation neu- 
trinos is sizable (sin^ 2^13 > 0.03) has further enhanced prospects for Hyper-Kamiokande (Hyper-K 
or HK) to discover CP asymmetry. 

Since 1970's, Grand Unified Theories (GUT or GUTs), which unify the strong and electroweak 
interactions and describe them arising from larger gauge symmetries like SU{5), have been exten- 
sively developed. Because leptons and quarks are often placed in the same multiplets, most GUTs 
allow baryon number violating interactions [THS]- Baryon number violating nucleon decays would 
constitute an extremely sensitive probe of BSM physics and the search for such a signal remains 
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one of the major endeavors in high energy physics. So far nucleon decays have escaped detection 
even by the world's largest nucleon decay detector, Super-K p^Hl3] . and so it must now be pursued 
by its successor, Hyper-K. 



1. Neutrino oscillations and CP violation 



The neutrino mixing matrix U |14l [T5] - often called the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) or Maki-Nakagawa-Sakata (MNS) mixing matrix - which translates neutrino mass eigen- 
states into flavor eigenstates (fe, f/^, z^r)^ = U{ui, 1/2, z^s)"^, is parameterized by three mixing angles 
(^12, ^23; ^13)) one Dirac CP phase 6, and two Majorana CP phases {021,0^1) in the three flavor 
neutrino framework [5] as 



U 



^ Uel Ue2 Ue3^ 

U^i Ufj_2 Uf^is 

\ Url Ur2 Ur3 J 

/l \ / 

C23 S23 

-S23 C23 J 



( Ci2 S12 ^ 



V 



-S12 C12 
1 



ll 



1 










,•"21 







e 2 











,■"31 

e 2 



1) 



ci3 sise * 
1 
y -siae*"^ ci3 J 

where Cjj and Sij represent smOij and cosOij, respectively. If massive neutrinos are Dirac particles, 
the Majorana phases are absorbed and only the Dirac phase can be responsible for CP violation 
in the lepton sector. Neutrino oscillation frequencies are determined by mass parameters, Am^]^ = 
1712 — m\ and Am32 = m\ — rriQ^, where mi,m2, and are the three mass eigenvalues. Because 
oscillation probabilities do not depend on the Majorana phases, the number of relevant oscillation 
parameters is six in total and the CP violation effect can be seen only via the Dirac phase 8 in 
neutrino oscillation experiments. Among these six parameters, 9i2 and Am-2x have been measured 
by solar and reactor neutrino experiments p^HT9] . ^23 and |Am32| (only its absolute value) have 
been measured by atmospheric neutrino |20| [2T] and accelerator neutrino experiments \22\ [23] . 

For many years, much effort has been devoted to measuring the last unknown mixing angle ^13 
by long baseline experiments [24j , atmospheric neutrino studies ^25j , as well as solar [ISj and reactor 
experiments p6]. Among them, the world's most stringent upper limit has been sin^ 2^13 < 0.15, 
which was obtained by the short baseline reactor neutrino experiment CHOOZ [27J. The quest has 
now been advanced by the T2K experiment announcing in June 2011 that a positive indication 
of — )• Uf. oscillations has been observed with the result that zero ^13 is disfavored at the 2.5a" 
level [1] . The significance becomes more than 3c7 |28j when a global fit is performed using available 
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data from MINOS, Super-K atmospheric u, solar, and reactor experiments in combination with 
the T2K observation. This indication - or evidence as it is called by the authors of j28j - would 
be a critical milestone in the overall neutrino physics program because a nonzero ^13 enables us to 
explore the 5 parameter. The magnitude of CP violation in i/i — )• f// and v'l — )• z/;/ oscillations is 
determined by the rephasing invariant: 

1 



Jcp = - cos sin 2^12 sin 2^23 sin 2^13 sin 6. 
8 



(2) 



Thus, the size of CP violation effect depends on the size of both ^13 and 6 values. 

The value of the CP phase is a question of general interest because (1) the 6 is last remaining 
oscillation parameter to be measured in the MNS matrix, and (2) a leptonic CP asymmetry is a 
leading candidate for understanding the origin of the existing matter-dominated (over antimatter) 
universe. One of the goals of the Hyper-K experiment is to measure the 6 parameter and thereby 
discover leptonic CP violation by using a J-PARC muon neutrino beam. We describe the expected 
sensitivity in Sec. Ill A 



(a) Normal hierarchy 



(b) Inverted hierarchy 



m 



m, 



m 



m, 



Am' ~ +2.4x10' eV- 




Am' ~+7.6HlO~'eV^ V, 



~ -2.4x10 -eV' 



FIG. 2. Two mass hierarchy cases. 

The mass hierarchy or the sign of Am32 is not yet known: whether 1713 is the heaviest mass 
eigenvalue (normal hierarchy, Am32 > 0) or is the lightest one (inverted hierarchy, Amgg < 0) 
remains to be experimentally determined in the future. Two hierarchy cases are illustrated in 
Fig. [2} Such an experimental determination would help in understanding the origin of neutrino 
masses, an explanation of which is expected to relate to BSM physics. Determination of the 
mass hierarchy is also important in terms of neutrino-less double beta decay rates. In the case 
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of inverted hierarchy, the effective neutrino mass m^^ is expected to be larger than 10 meV; 
this could be reachable by various proposed detectors designed to test the Majorana nature of 
neutrinos |29j . The suggested sin^ ^13 of a few xlO^^ would also open the possibility of mass 
hierarchy determination by measuring atmospheric z^e appearance via a resonance effect arising 
from passing through the Earth's matter. High statistics atmospheric neutrino data would also 



provide additional information on 5 and the ^23 octant, as discussed in detail in Sec. IIIB If we are 
fortunate to observe the neutrino burst from a supernova explosion, the mass hierarchy could be 



determined by observing time variation of the neutrino spectrum (see Sec. HIE 1 for the details). 



2. Nucleon decays 



Motivations for unifying the strong, weak, and electromagnetic forces are very strong. Grand 
Unifies Theories naturally explain the experimental observation that electrical charges are quan- 
tized as Qproton + Qe- = with better than 10~^^ accuracy. GUTs provide a simple organization 
of quarks and leptons by which the quantum numbers of quarks and leptons may be naturally un- 
derstood. GUTs are also supported by the noted unification of three gauge coupling constants at 
the very high energy scale of Mqut ~ 10^^ GeV. They have further support from the high energy 
physics scale implied by the tiny neutrino mass scale observed in neutrino oscillation experiments. 

In the SM, protons are practically stable particles because the decay rate due to the chiral 
anomaly is unobservably small [30]. However, GUTs generally yield a prediction of spontaneous 
proton or bound nucleon decay with experimentally testable decay rates. Therefore nucleon decay 
experiments are well motivated to test various proposed unified models, detect evidence of a new 



{a) p^e^ +Jt^ 



u 



P< u 

d 



X 



(b) p^y + K^ 



d 
d 




FIG. 3. Examples of Feynman diagrams for particular nucleon decay modes. The left panel shows a proton 
decaying into a positron and a neutral pion via exchange of a super-heavy gauge boson X. The right one 
shows the SUSY-favored decay mode, proton decay into a neutrino and a mediated by a colored Higgs 
and other SUSY particles. 
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TABLE III. Proton lifetimes predicted by various GUTs and SUSY-GUTs. 



Model 


Decay Mode 


Lifetime Prediction (years) 


Minimal SU{5) 


P- 




1028.5 , 


-1031-5 [H] 


Minimal 50(10) 


P- 


^e+7rO 


1030- 


1040 1311 


Minimal SUSY SU{5) 


P- 


> 1^K+ 


< 1030 


[33] 


SUGRA SU{5) 


P- 


> VK+ 


1032- 


1034 34J 


SUSY 50(10) 


P- 


> VK+ 


1032- 


1034 


SUSY 50(10) 


P- 


^e+7rO 


< 5.3 X 1034 36^ 



paradigm, and then to pin down the details of the successful theory by measuring the branching 
ratio of each nucleon decay mode. 

The favored decay mode in prototypical GUTs based on SU (5) symmetry is p — t- e"*" + vr'^ 
mediated by super-heavy gauge bosons as shown in Fig. [3} On the other hand, GUTs incorporating 
supersymmetry (SUSY-GUTs) suppress the decay mode p — >• e"^ + vr'^ while favoring p — )• z7 + 
via dimension five operator interactions involving exchange of a heavy color triplet Higgsino as 
shown in the right panel of Fig. [3j Proton lifetimes predicted by some GUTs and SUSY-GUTs 
are listed in Table Iml 

For three decades now, the experimental pursuit of nucleon decay has been led by water 
Cherenkov detectors. For example, increasingly stringent limits on the p — )• e"*" + vr'' decay mode 
have been set by a series of experiments employing this technology: 2.6 x 10^^ years by Kamiokande 
in 1989 [33, 8.5 X 10^^ yg^rs by 1MB in 1999 I38j, and 8.2 x 10^^ years by Super-Kamiokande in 
2009 [13]. Due to these null observations, minimal SU{5) is now an excluded model. In addition, 
minimal SUSY SU{5) is also considered to be ruled out |39| by the results of p — )• VK^ searches, for 
which water Cherenkov detectors also provide best limits p!2] . It is certain that a water Cherenkov 
detector with 25 times the fiducial volume of Super-Kamiokande will provide further stringent 
tests of various GUTs. Such a large detector will have a good chance to discover nucleon decay 
phenomena and would then have the potential to provide important constraints on the emergent 
unified theory by measuring the relative rates of a variety of nucleon decay modes. The expected 



sensitivities for various nucleon decay modes are described in Sec. Ill C 
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3. Astrophysical objects 

Hyper-Kamiokande is also capable of observing neutrinos from various astrophysical objects. 
One of advantages of the detector is that its energy threshold is as low as several MeV; this enables 
us to reconstruct neutrinos from the Sun and supernova on an event-by-event basis. 

For 40 years, many efforts have been made to computationally simulate core collapse supernova 
explosions, but the simulated bounce shock waves often stalled or required special assumptions to 
be made. It would seem that the details of the supernova explosion mechanism are still lacking. 
For instance, one important effect is thought to be the matter reheating process from neutrinos 
trapped in the iron core. High statistics observations of neutrinos from a supernova would provide 
precious inside information of the explosion mechanism. If a core collapse supernova explosion 
were to take place near the center of our Galaxy, Hyper-K would observe as many as 250,000 
neutrino interactions. There is also a good chance to explore the mechanism of neutron star or 
black hole formation through the fine-grained, sub-millisecond time structure details of such a 
burst. Meanwhile, while waiting for a nearby explosion to occur, the continuous flux of supernova 
neutrinos from all past core collapse explosions will guarantee a steady accumulation of valuable 
astrophysical data. 

Thanks to its good low energy performance for upward-going muons, Hyper-K has a larger effec- 
tive area for upward-going muons below 30 GeV than do cubic kilometer-scale neutrino telescopes 



(see Fig. 65 in Sec. HIE). Additionally, fully contained events in Hyper-K have energy, direction, 
and flavor reconstruction and resolutions as good as those in Super-K. This high performance will 
be useful for further background suppression or studies of source properties. For example, the 
detector is extremely sensitive to the energy range of neutrinos from annihilations of light (below 
100 GeV) WIMP dark matter, a region which is suggested by recent direct dark matter search 
experiments. 



Sensitivities and expected physics outputs from astrophysics studies are discussed in Sec. |IIID 
andlHEl 



B. Water Cherenkov technique 

The technique of ring imaging water Cherenkov detection, which the Hyper-Kamiokande detec- 
tor is based on, has been successfully employed by several neutrino and nucleon decay experiments 
over the last few decades. In Japan, the Kamiokande detector (1983-1996, 3 kiloton) and Super- 
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FIG. 4. Event displays of observed data from Super-Kamiokande's phase-I period. Shown are unrolled views 
of a single-ring electron- like event (top left), a single-ring muon-like event (top right), and a two- ring event 
(bottom). Colored points indicate the number of detected photo-electrons in each photomultiplier tube. 

Kamiokande (in operation since 1996, 50 kiloton) have achieved several important scientific results, 
notably detection of extraterrestrial neutrinos from the Sun [30] and Supernova 1987a |4HI42j. and 
discovery of neutrino flavor mixing and neutrino mass [6l |l3] . In the K2K long baseline neutrino 
oscillation experiment, Super-K and a one kiloton water Cherenkov detector (IKT) provided indis- 
pensable data on the neutrino beam flux and its energy spectrum at the neutrino production site 
(using IKT) and a location 250 km farther away (using Super-K) [22]. Super-K again is playing 
the role of the far detector in the ongoing T2K experiment which reported an indication of — )• i/^ 
oscillations in June 2011 [1]. 

Relativistic charged particles, produced via neutrino interactions or possible nucleon decays. 
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TABLE IV. Expected detector performance of Hyper-Kamiokande. Further details can be found in Sec. Ill 



Resolution or Efficiency 



Vertex resolution @ 500MeV/c 
@ 5GeV/c 



28 cm (electron) / 23 cm (muon) 
27 cm (electron) / 32 cm (muon) 



Particle ID @ 500MeV/c 
@ 5GeV/c 



98.5 ± 0.6 % (electron) / 99.0 ± 0.2 % (muon) 
99.8 ± 0.2 % (electron) / I00lo° % (muon) 



Momentum resolution @ 500MeV/c 
@ 5GeV/c 



5.6 % (electron) / 3.6 % (muon) 
2.0 % (electron) / 1.6 % (muon) 



Electron tagging from 500MeV/c /x+ decays 
from 5GcV/c /i+ decays 



98 % 
58 % 



J-PARC I'e signal efficiency 



64 % (nominal) / 50 % (tight) 



J-PARC i^fj CC background rejection 



>99.9 % 



J-PARC ly 7r° background rejection 



95 % (nominal) / 97.6 % (tight) 



p + TT^ efficiency (w/ tt^ intra- nuclear scattering) 

atmospheric v background 



45 % 

1.6 events/Mton/year 



p v + efficiency by prompt 7 tagging method 
atmospheric v background 



7.1 % 

1.6 events/Mton/year 



p^ D + K+,K+ 



7r° efficiency 



atmospheric ly background 



6.7 % 

6.7 events/Mton/year 



Vertex resolution for 10 MeV electrons 



90 cm 



Angular resolution for 10 MeV electrons 



30° 



Energy resolution for 10 MeV electrons 



20 % 



radiate Cherenkov photons while passing through the water [Hj. These photons' production angle 
compared to the charged particle's direction is given by cos 0c = l/nj3, and 6c is 42 degrees in the 
case that the particle velocity j3 = v / c = 1 and the refractive index of the medium is n = 1.33 (i.e., 
water). By detecting the spatial and timing distributions of the Cherenkov photons with an array 
of single photon capable sensors, the parent charged particle's position, direction, and energy can 
be reconstructed. Ring imaging water Cherenkov detectors are also able to discriminate charged 
muons (/U^) from electrons, positrons, and gammas (e^ / 7), as the later group of particles induce 
electromagnetic showers in the medium which modify the resulting ring patterns. Figure |4] shows 
event displays of observed data from Super-Kamiokande's phase-I data-taking period (1996-2001). 
The fraction of mis- identification of muons (electrons) as electrons (muons) is as low as 1% as shown 
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in Table IV Moreover, the high efficiency of tagging muon decay electrons provides additional 
information on the particle type with which to further purify muon or electron samples. Thanks to 
their excellent particle identification capabilities, water Cherenkov detectors can determine neutrino 
flavor compositions in atmospheric and accelerator neutrinos and thereby detect neutrino flavor 
transitions or oscillations. Other critical features of water Cherenkov detectors are an excellent 
rejection efficiency for vr'' background interactions as well as a good signal efficiency for — t- 
oscillation studies using the J-PARC neutrino beam. These particle ID and tracking capabilities 
are also indispensable for identifying possible nucleon decay candidates while at the same time 
minimizing physics backgrounds to nucleon decay arising from atmospheric neutrino interactions. 



Table IV presents a summary of Hyper-Kamiokande's high-energy performance relevant for 
nucleon decay searches as well as atmospheric and accelerator neutrino oscillation experiments in 
which particle energies are above ~ 100 MeV. Detector performance for low-energy events such as 
solar neutrinos or supernova neutrinos are also shown. Details of data analyses and the physics 



sensitivities of the Hyper-K detector are discussed in Sec. Ill 
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II. EXPERIMENTAL SETUP OVERVIEW 

This section explains the basehne design of the Hyper-Kamiokande detector, i.e. candidate 
site, cavern and water tank, water purification system, photosensors, data acquisition system, and 
cahbration system. Necessary technology for the detector has been already established and the 
detector could be constructed and fully operational in 10 years. 

A. Site, caverns, and tanks 




FIG. 5. The candidate site map. The site is located 8 km south of the Super-K site as shown in the left 
panel. The map of the Tochibora mine is shown in the right panel. 

Hyper-Kamiokande is to be the third generation water Cherenkov detector in Kamioka designed 
for nucleon decay searches and neutrino studies. Its total volume of one million tons would be 
about 20 times larger than that of the world's largest currently operating underground detector, 
Super-Kamiokande. The Hyper-K detector candidate site, located 8 km south of Super-K, is in 
the Tochibora mine of the Kamioka Mining and Smelting Company, near Kamioka town in Gifu 
prefecture, Japan, as shown in Fig. [5j The experiment site is accessible via a drive-in, ~2.6 km long, 
(nominally) horizontal mine tunnel. The detector will lie under the peak of Nijuugo-yama, having 
648 meters of rock or 1,750 meters- water-equivalent (m.w.e.) overburden, at geographic coordinates 
Lat. 36°21'08.928"N, Long. 137°18'49.688"E, and an altitude of 508 m above sea level. The rock 
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FIG. 6. Cross section view of the Hyper-Kamiokande detector. 



wall in the existing tunnels and sampled bore-hall cores are dominated by Hornblende Biotite 
Gneiss and Migmatite in the state of sound, intact rock mass. This is desirable for constructing 
such unprecedented large underground cavities. The site has a neighboring mountain, Maru-yama, 
just 2.3 km away, whose collapsed peak enables us to easily dispose of more than one million m^ 
of waste rock from the detector cavern excavation. The site also has the benefit - well-suited for 
a water Cherenkov experiment - of abundant, naturally clean water located nearby. More than 
13,000 mVday (i.e., one million tons per ~80 days) will be available. The Mozumi mine where 
the Super-K detector is located is another candidate site for which less geological information is 
available at this moment. 

In the base design, the Hyper-K detector is composed of two separated caverns as shown in 
Fig. [T| each having a egg-shape cross section 48 meters wide, 54 meters tall, and 250 meters long 
as shown in Fig. [6j [7j and [8j These welded polyethylene tanks are filled up to a depth of 48 m 
with ultra-pure water: the total water mass equals 0.99 million tons. The 5 mm thick polyethylene 
sheet which forms the water barrier has been used in waste disposal sites, and based upon that 
experience its expected lifetime is longer than 30 years. We expect a much longer lifetime than 
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FIG. 7. Cross section view of the outer detector. 



30 years, as the underground experimental site is free from the sunshine which degrades plastics 
on the Earth's surface. This polyethylene lining sheet will be placed upon a concrete base to 
achieve a strong, water tight structure as is shown in Fig. [7) Detailed design of the structure 
of the concrete layer, connection with the lining sheet, and the PMT support structure is yet to 
be made. Each cavern will be optically separated by segmentation walls located every 49.5 m 
to form 5 (in total 10) compartments as shown in Fig. [9| such that event triggering and event 
reconstruction can be performed in each compartment separately and independently. Because the 
compartment dimension of 50 m is comparable with that of Super-K (36 m) and is shorter than the 
typical light attenuation length in water achieved by the Super-K water filtration system, (> 100 
m @ 400 nm), we expect that detector performance of Hyper-K will be basically the same as 
that of Super-K. The water in each compartment is further optically separated into three regions. 
The inner region has a barrel shape 42 m in height and in width, and 48.5 m in length, and 
is viewed by inward- facing array of 20-inch diameter photomultiplier tubes (PMTs). The entire 
array consists of 99,000 HAMAMATSU R3600 hemispherical PMTs, uniformly surrounding the 
region and giving a photocathode coverage of 20%. The PMT type, size, and number density are 
subject to optimization. An outer region completely surrounds the 5 (in total 10) inner regions 
and is equipped with 25,000 8-inch diameter PMTs. This region is 2 m thick at top, bottom, and 
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FIG. 8. Profile of the Hyper-K detector. The left panel shows PMT arrays and the support structure for 
the inner and outer detectors. The right panel shows segmentation of the detector. Each quasi-cylindrical 
tank lying horizontally is segmented by intermediate walls into five sub-detectors. 
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barrel sides, except at both ends of each cavern, where the outer region is larger than 2 m due 
to rock engineering considerations. A primary function of the outer detector is to reject entering 
cosmic-ray muon backgrounds and to help in identifying nucleon decays and neutrino interactions 
occurring in the inner detector. The middle region or dead space is an uninstrumented, 0.9 m 
thick shell between the inner and outer detector volumes and the stainless steel PMT support 
structure is located in this region. Borders of both inner and the outer regions are lined with 
opaque sheets. This dead space, along with the outer region, acts as a shield against radioactivity 
from the surrounding rock. The total volume of the inner region is 0.74 million tons and the total 
fiducial volume is 10 times 0.056 = 0.56 million tons. The fiducial volume is defined as the region 
formed by a virtual boundary located 2 m away from the inner PMT plane (Fig. ^ right panel). 



Figure 10 shows a schematic view of the cavern construction. After constructing access tunnels, 
approach tunnels, and belt conveyor tunnels, construction of two water tank caverns will start 
with the excavation of arches, and then excavation will proceed from top to bottom for each layer 
utilizing the bench cut method. The majority of waste rock is planned to be transported to the 
top of the adjacent mountain Maru-yama via conveyors in the declined conveyor tunnels. 

The estimated cosmic-ray muon rate around the Hyper-K detector candidate site is shown 



in Fig. 11, The range of estimated muon fluxes at the detector position is 1.0 ~ 2.3 x 10^ 
sec~^cm~^ which is roughly 10 times higher than the flux at Super-K's location (0.13 ~ 0.14 x 
10^^ sec~^cm~^). The expected deadtime due to these muons is less than 1% and negligible for 



24 II EXPERIMENTAL SETUP OVERVIEW 




FIG. 10. Schematic view of the cavern construction. 



nucleoli decay searches, atmospheric neutrino studies, and long baseline experiments. However, 
the radioactive spallation products produced by the higher cosmic-ray muon rate would increase 
background levels for low-energy physics targets such as solar or supernova neutrinos. An estima- 
tion of the possible background increase by the higher cosmic-ray muon rate - and its impact - 
will be discussed in the Sec. IIIIDl 



B. Water purification system 

Needless to say, water is the target material and signal-sensitive medium of the detector, and 
thus its quality directly affects the sensitivity. In order to realize such a huge Cherenkov detector, 
achieving good water transparency is the highest priority. In addition, as radon emanating from 
the photosensors and detector structure materials is the main background source for low energy 
neutrino studies, an efficient radon removal system is indispensable. 

In Super-Kamiokande the water purification system has been continually modified and improved 
over the course of SK-I to SK-IV. As a result, the transparency is now kept above 100 m and is very 
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FIG. 11. Muon flux calculated by MC simulation around the Hyper-K candidate location at Lat. 
36°21'08.928"N, Long. 137°18'49.688"E, and altitude 508 m above sea level. Horizontal and vertical axes 
represent distance from the candidate location in m in the North-South and West-East direction, respec- 
tively. Text and color represents expected muon flux in units of 10~^sec~^cm~^. Two boxes represent the 
candidate locations of the two Hyper-K detector modules. The calculation has been performed based on 
the predicted muon flux at the Earth's surface, attenuation by the rock around the detector. The depth in 
the mountain is calculated using a 50 m mesh tomographic map. A rock density of 2.7g/cm^ is assumed in 
the attenuation calculation. 

stable, and the radon concentration in the tank is held below 1 mBq/m^. Following this success, 
the Hyper-Kamiokande water system design will be based on the current Super-Kamiokande water 
system. 

Naturally, ever-faster water circulation is generally more effective when trying to keep huge 
amounts of water clean and clear, but increasing costs limit this straightforward approach so a 
compromise between transparency and recirculation rate must be found. In Super-Kamiokande, 
50,000 tons of water is processed at the rate of 60 tons/hour in order to keep the water transparency 
(the attenuation length for 400 nm-500 nm photons) above 100 m, and 20 m'^/hour of radon free air 
is generated for use as a purge gas in degas modules, and as gas blankets for both buffer tanks and 
the Super-Kamiokande tank itself [l5]. For the 0.99 million tons of water in Hyper-Kamiokande, 
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FIG. 12. Current air & water flow design of the Hyper-Kamiokande water system. 



these process speeds will need to be scaled-up to 1200 m'^/liour for water circulation and 400 
m^/hour for radon free air generation. 



Figure 12 shows the current design of the Hyper-Kamiokande air and water purification system. 
With these systems, the water quality in Hyper-Kamiokande is expected to be same as that in 
Super-Kamiokande. 



B Water purification system 
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C. Photosensors 

In addition to clean water, photosensors are the key component of a water Cherenkov detector. 
The intrinsic characteristics of such a detector impose a number of requirements on the crucial 
photodetectors. 

Of course, a water Cherenkov detector uses the Cherenkov photons to identify and to reconstruct 
the directions, energies, and locations of the particles which produce them. The total number of 
Cherenkov photons produced has a linear correlation with the momentum of the particle. Therefore, 
the total number of Cherenkov photons of a low energy particle - for example, a recoil electron from 
a solar neutrino elastic scattering interaction - are few, and so the photosensors in the detector 
are required to have high efficiency and good resolution at the level of a single photon. On the 
other hand, several hundred photons hit each PMT when high energy particles - for example, the 
particles generated by an inelastic atmospheric neutrino interaction - pass through the detector. 
Therefore, the photosensors are also required to have wide dynamic range and good linearity. 
Otherwise, it is not possible to achieve good particle energy resolution. 

The location of the interaction vertex is reconstructed using photon arrival timing information 
at each PMT. Therefore, good timing resolution of the photosensors is essential, and the jitter of 
the transit time is required to be small, i.e. less than 3 nanoseconds for a single photon. 



1. Photosensors in the baseline design 

In the baseline design, we have selected the 20-inch PMT ( R3600 ) used successfully in Super- 
Kamiokande as the primary sensor candidate; this type of PMT is already known to satisfy the 
requirements above. Moreover, it been operated for more than 15 years in Super-Kamiokande and 
thus not only the performance characteristics but also the long-term stability of this type of sensor 
are well understood. The specifications of the 20-inch PMT is summarized in Table [V] while the 



single photoelectron distribution for this PMT is shown in Fig. 13 

Based on studies of the expected physics performance of Hyper-Kamiokande, it is necessary to 
cover 20% of the detector wall with these PMTs. This means that in order to meet this physics- 
based requirement the total number of PMTs needed will be about 99,000. 

In order to use R3600 in the Hyper-K detector, it is necessary to equip each PMT with a 
protective case. Otherwise, a collapse of one of the sensors could generate a shock wave large enough 
to break the surrounding PMTs, potentially leading to a chain reaction resulting in catastrophic 
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Shape 


Hemispherical 


Photocathode area 


50 cm diameter 


Window material 


Pyrex glass (4^5 mm) 


Photocathode material 


Bialkah (Sb-K-Cs) 


Quantum efhciency 


22% at A = 390 nm 


Dynodes 


11 stage Venetian blind type 


Gain 


10^ at - 2000 V 


Dark current 


200 nA at 10^ gain 


Dark pulse rate 


3 kHz at 10^ gain 


Cathode non-uniformity < 10 % 


Anode non-uniformity 


< 40% 


Transit time 


90nsec at 10 gain 


Transit time spread 


2.2nsec {la) for 1 p.e. equivalent signals 


Weight 


13 kg 


Pressure tolerance 


6 kg/cm^ water proof 


TABLE V. 


Specifications of 20-inch PMT. 
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FIG. 13. Single photoelectron distribution of a typical 20-inch PMT. 



damage to the detector. In order to avoid this kind of severe accident, each PMT will be equipped 
with a protective, shock-absorbing case specifically designed to prevent a collapse of a single PMT 
from emitting such a fratricidal shock wave. 

Currently, the PMTs in Super-K are all contained within 10 mm thick FRP cases in back and 



the 13 mm thick acrylic covers in front as shown in Fig. 14 This type of cover is known to satisfy 
the minimum requirements of the Hyper-K detector. 



C Photosensors 
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2. Possible alternative PMTs and other R&D items 

If photosensors with higher efficiency are available, it would be possible to reduce the number 
or the size of the photosensors in Hyper-Kamiokande. There are several ongoing efforts to build 
high efficiency sensors. One approach is a small diameter PMT ( ~ 10 inch ) with high quantum 
efficiency (QE). Its QE is expected to be ~ 40%, which is almost double that of the R3600. Another 
path would be to utilize a new type of sensor, called a Hybrid PMT. The Hybrid PMT is a sensor 
which utilize an avalanche photo diode (APD) instead of dynodes; it also has higher single photon 
sensitivity and much higher timing resolution. 

Of course, it is necessary to increase the number of PMTs by 1.4 if 10 inch PMTs with doubled 
QE are used. However, the relative smallness of the sensor will be a great help in safe handling 
during the construction. Also, protective cases would presumably be less expensive compared to 

the ones for R3600 for two reasons: smaller size means less case material, and a smaller PMT 
means considerably less stored energy to contain in case of breakage, so the cases could likely be 
thinner as well as smaller. 

There have also been several attempts to improve the photon collection efficiency with special 
lens systems or mirrors attached to the existing sensors. It is well to study this kind of idea, but 
whether or not it can be made to work with this type of detector remains to be seen. Often such 
optical approaches to increase light collection come with a price; for example, the loss of effective 
fiducial volume due to limited angles of acceptance, or decreased light arrival timing resolution. 

The other important R&D item is the design of the protective case. The current design used 
in Super-K is known to meet the requirements, but as a retrofitted design it was forced to satisfy 
several additional requirements only applicable to Super-Kamiokande. In Super-K it was not 
possible to change the mounting method of the PMTs, and the quantity of cases produced during 
each of two production runs was just 5,000 or so, not enough to benefit from economies of scale 
or automated processes. In the case of Hyper-Kamiokande it will be possible to design the frame 
mounts of the photosensors and the protective cases simultaneously, and it will be possible to 
simplify the design of the case as well as the mounting methods. As for the choice of material used 
in the protective cases in Super-K, there were restrictions due to the limited quantity required. 
Based on the Super-K experience, FRP is now known to contain non-negligible amounts radioactive 
substances. In Hyper-K we will produce at least 20 times as many cases, and thus it should not be 
necessary to stick to FRP but instead to use a low-activity material like stainless steel or acrylic. 
The design of the PMT frame mounts and protective cases has been started in coordination with 
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FIG. 14. PMT case attached a 20-inch PMT. These were added to all Super-K phototubes in 2002. 



Charge dynamic range ^0.2 pC ^ 2500 pC 

Integral charge non-linearity less than 1% 

Charge resolution 0.2 pC 

Discriminator threshold (noise level) less than ImV 

Timing resolution much better than Ins @ 2pC 

Power consumption per channel less than IW/ch 

Data through rate lOOkHz for 10 seconds 

lOkHz as nominal dark rate 

Channel deadtime Smaller than ~ 500ns 

TABLE VI. Requirements for the front-end electronics. 



the design of the detector structure. 



D. DAQ system 



The required specifications for the front-end electronics are basically the same as the current 
Super-K modules, assuming that the Hyper-K PMT gain and the noise levels are also similar. 



These specs are as summarized in Table VI 



In Super-K we have developed a free-running, triggerless (at the hardware level) DAQ system. 



The schematic diagram of this DAQ is shown in Fig. 15 

The front-end boards digitize any signals which exceed the discriminator threshold and send out 
the charge and timing information to the front-end computers. The front-end computers collect hit 
information from 30 electronics modules each, corresponding to ~ 500 PMTs, sort the hit data in 
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To offline system 
( outside of ttie mine ) 

FIG. 15. Schematic diagram of the data acquisition system. 

the order of timing and send it out to the merger computers. The merger computers each collect all 
the hit information from a compartment and apply a software trigger to remove noise hits. Then, 
the organizer computer collects the event from the mergers, eliminates overlapping events, and 
send them to the offline computer system. This system reads out all the digitized hit information 
from the PMTs and selects the events with software. This system has been working without any 
problem in Super-K for more than 2 years. Therefore, it is safe to say that there are no serious 
technical difficulties in preparing the DAQ system for this new detector with the currently available 
and well established tools at hand. The expected data rate from the entire Hyper-K detector is 
about 5 GB/s before applying the software trigger. After the first level software trigger, it will 
be reduced by 1/50 and another factor of five reduction could be achieved with Super- K-style 
intelligent vertex fitters applied in the offline computer system. In the end, the data rate written 
to disk is expected to be less than 100 MB/sec in total. 

1. R&D items for the DAQ system 

A possible difficulty which could arise would involve physically running over 100,000 cables from 
the PMTs - arranged as they must be throughout the immense tank - into the DAQ system. The 
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level of challenge will heavily depend on both the method selected to mount the sensors to their 
supporting framework as well as how the cables are connected to the sensors, but it will surely not 
be a simple task and will contribute a significant amount of time to the total detector construction 
schedule. The great mass of cables inside the detector will also inevitably create some insensitive 

regions, particularly around the cable feed through leading out of the tank. Moreover, the length 
of the cables could be longer than 100 m, some 50% longer than in Super-K. This may degrade the 
quality of the signal. 

Therefore, it is worthwhile to consider putting the front-end electronics inside the detector itself. 
If it is possible to place the electronics modules close to the photosensors in the water then not 
only will the cable lengths be shorter, but also the number of cables needing to pass through the 
detector surface will be reduced. We will need to consider the resulting temperature rise if we 
put electronics modules in the water. However, the current power consumption of the digitization 
module and the high voltage supplies is a few watts per channel; this should be acceptable even if 
we use more than 100k sensors. 

There are still several R&D items necessary to allow placing the electronics modules in the 
water. These include the fault-tolerance in the circuit, a water-proof case and its connectors or 
cable feed through, and a proper design of the heat exchanger from the circuit. The high voltage 
power supply for the photosensor will also necessarily be placed in the same housing. There already 
exist possible solutions for each such item, and thus what we have to do from an R&D standpoint 
is evaluate each implementation and select the proper ones considering both quality and cost. 

E. Detector calibration 

In Super-Kamiokande, various kind of calibrations have been performed. A nitrogen+dye laser 
- which emits pulses of light with a wavelength of 396 nanometers, duration of < 4 nanoseconds, 
and energy of 300 microjoules - has been used to calibrate the timing response of PMTs as well 
as measure the charge linearity of PMTs + DAQ electronics. The laser light is feed to a diffuser 
ball in the Super-K tank through an optical fiber with a large (200 ;um) diameter core. 

A nitrogen laser with 337 nm wavelength and laser diodes with 375 nm, 405 nm, 445 nm, 
and 473 nm wavelengths have been used to measure water attenuation and scattering parameters. 
Through optical fibers, these five laser light sources power eight light injectors which are mounted 
in various positions inside the Super-K tank. A computer-controlled optical switch directs the laser 
light to the various injectors, allowing the measurement of absorption and scattering parameters 
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for light in the water with wavelength dependence taken into account. 

A xenon lamp which has stable( ~5% ) light intensity has been used to calibrate PMT gain 
and to monitor the position dependence of water quality. UV light from the xenon lamp travels 
through an optical fiber with very large (1000 iim) diameter core, and into a scintillator ball sitting 
inside the Supcr-K tank. 

An surplus medical grade electron linear accelerator (LINAC) system is permanently mounted 
near the top of the Super-K tank. It allows single 5^16 MeV electrons to be injected into various 
positions in the Super-K tank through an evacuated beam pipe and bending magnets which are 
positioned as needed to guide the beam. At one electron every second this is likely the lowest flux 
particle accelerator in the world, but it has proven to be our most accurate low energy calibration 
source. Through the use of a germanium counter the beam energy is known to within 20 keV. 

A deuterium-tritium (DT) fusion generator which produces gamma and electrons via + 

"^He + n, ^*^0(n,p)^^N, ^^^N 7(6.1MeV), e(4.3MeV) is able to be lowered into the Super-K 
tank when a gamma or electron source in this energy range is needed. As it is much less labor 
intensive to operate than the LINAC, the DT generator is run about four times as frequently - a 
few times each year - to ensure that the detector response is stable and its response uniform to 
low energy events. 

A gamma ray source powered by neutrons emitted by from ^^^Cf produces ~9 MeV gamma 
rays from nickel via Ni(n, 7)Ni. Like the DT generator, it is able to be lowered directly into the 
Super-K tank. All of these low energy gamma ray and electron sources are used to calibrate vertex 
and energy determination, primarily for use in solar and supernova neutrino analysis. 

Cosmic ray muons which either pass all the way through or stop somewhere in the tank, as well 
as the decay electrons produced by these stopping muons are used to calibrate energy reconstruction 
for more energetic events. The cosmic ray muons are also used as a natural source of Cherenkov 
light for monitoring water transparency. 

The laser and xenon light sources (and of course the cosmic ray muons) are used even during 
normal data taking by keeping the diffuscr ball, scintillator ball, and light injectors in the tank 
and scheduling light emission times of each light source at a low but steady rate(~l Hz in total) 
to minimize dead time of the observation. 

The LINAC and radioactive calibration sources cannot be used during normal data taking. 
Therefore, for calibration work with these sources; for high rate data taking of laser and xenon; 
and to install, uninstall, or move calibration sources in the tank, we need to stop normal data 
taking, leading to observational dead time. This is purposely kept to a minimum, primarily to 
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reduce the risk of missing a once-in-lifetime burst of neutrinos from a galactic supernova, but also 
to avoid wasting expensive accelerator beam time whenever a long baseline experiment is running. 
Calibrations can often be planned around scheduled accelerator downtime, but the risk of missing 
a supernova is obviously always present when the detector is turned off. 

Almost all of these calibration sources will prove indispensable for Hypcr-K as well to maintain 
the best performance of the new detector. Since in the current baseline design of Hyper-K the 
two tanks are separated into 10 individual compartments, it would require - roughly speaking - 
10 times more time and manpower to perform the same calibration work. In addition, due to its 
egg-shaped cross section, the dome above the Hyper-K tank will be quite a bit narrower compared 
to the maximum tank width. This may make calibrations near the wall of PMTs challenging, as 
will the bowed walls of the detector. 

Therefore, calibration methods should be carefully discussed during the detector design, and 
the detector should have dedicated, automated systems for placing various calibration sources at 
desired positions within the tank. Due to the greater size and more complex geometry of Hyper-K, 
the positioning of calibration sources may need to be controlled by more advanced methods than 
the vertical wire drop lines used in Super-K. Convenient methods like fixed rails, wire guides, or 
even submarine robots installed during detector construction should be considered. As for the 
large, immovable instruments like the LIN AC, we may need to consider installing a small network 
of simple beamlines to distribute particles to each compartment. 
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III. PHYSICS POTENTIAL 

A. Accelerator based neutrinos 

1. CP asymmetry measurement in a long baseline experiment 

If a finite value of 613 is discovered by the ongoing and near- future accelerator and/or reactor 
neutrino experiments |46H50j . the next crucial step in neutrino physics will be the search for CP 
asymmetry in the lepton sector. A comparison of muon-type to electron- type transition probabil- 
ities between neutrinos and anti-neutrinos is one of the most promising methods to observe the 
lepton CP asymmetry. Recent indication of a nonzero, rather large value of ^13 [1] makes this 
exciting possibility more realistic with near-future experiments such as Hyper-Kamiokande. 

In the framework of the standard three flavor mixing, the oscillation probability is written using 



the parameters of the MNS matrix (see Sec. I A 1 ), to the first order of the matter effect, as jSlj : 



^ ^e) = ^Cf^SfA • sin^ A31 

-h8Cf3S'i25'i3S'23(Ci2C23 COS 6 - S12S13S23) ' COS A32 • sin A31 • sinA2i 
-8Ci3Ci2C235'i25'i3523sin5 • sin A32 • sin A31 • sinA2i 
+4.Sl2Cf3iCf2Cis + Sf2Si3Sl^ - 2C12C23S12S23S13 cosS) ■ sin2 A21 
-8CI3SI3SI3 ■ ^(1 - 2SI3) ■ cos A32 • sin A31 

+8C!3S!3Sl3-^{l - 2S!3) ■ sin' A31, (3) 

where Cij, Sij, Ajj are cosOij, sinOij, Amfj L/AE,^, respectively, and a[eV^] = 7.56 x 10^^ x 
p[g/cm'^] X Ejy[GeY]. The parameter 6 is the complex phase that violates CP symmetry. The 
corresponding probability for z7^ — )• transition is obtained by replacing 5^—6 and a — )• —a. 
The third term, containing sin (5, is the CP violating term which flips the sign between 1/ and P 
and thus introduces CP asymmetry if sin 5 is non-zero. The last two terms are due to the matter 
effect; caused by coherent forward scattering in matter, they produce a fake (i.e., not CP-related) 
asymmetry between neutrinos and anti-neutrinos. As seen from the definition of a, the amount 
of asymmetry due to the matter effect is proportional to the neutrino energy at a fixed value of 
L/E,. 



Figure 16 shows the — )• Ue and z?^ — )• i^e oscillation probabilities as a function of the true 
neutrino energy for a baseline of 295 km. The parameters other than ^13 and 6 assumed in 
this section are summarized in Table VII The value of sin^ ^23 is set to the maximal mixing, as 
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neutrino anti-neutrino 
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FIG. 16. Oscillation probabilities as a function of the neutrino energy for — >■ v^. (left) and — >■ (right) 
transitions with L=295 km and sin^ 2^13 = 0.1. Black, red, green, and blue lines correspond to ^ = 0, W, tt. 



and — ^TT, respectively. Other parameters are listed in Table VII Solid (dashed) line represents the case for 
a normal (inverted) mass hierarchy. 




FIG. 17. Oscillation probability of i/^ ^ i^e as a function of the neutrino energy with a baseline of 295 km. 
Left: sin^ 2013 — 0.1, right: sin^ 2(?i3 = 0.01. 5 ^ and normal hierarchy is assumed. Contribution from 
each term of the oscillation probability formula is shown separately. 



TABLE VII. Parameters other than 6*13 and S assumed in this section. 



Name 


Value 


L 


295 km 




T.exlO"'^ 


1 Arrival 


2.4x10^3 eV^ 


sin^ 6*12 


0.31 


sin^ 6*23 


0.5 


Density of the earth (p) 


2.6 g/cm'"^ 
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suggested by the current best fit values of atmospheric and accelerator-based experiments [201 - 
23| [52] . The density of earth is based on [53] . Other parameters are based on the global fit |28] . 
The cases for 5 = 0, ^it,tt, and — ^vr, are overlaid in Fig. 



16 



Also shown are the case of normal 



mass hierarchy (Amgg > 0) with solid lines and inverted mass hierarchy < 0) with dashed 

lines. The oscillation probabilities depend on the value of 5, and by comparing the neutrinos and 
anti-neutrinos, one can see the effect of CP violation. There are sets of different mass hierarchy 
and values of 6 which give similar oscillation probabilities. This is known as the degeneracy due 
to unknown mass hierarchy and may introduce a fake solution if we do not know the true mass 
hierarchy. The effect of this degeneracy on the measurement of CP asymmetry will be discussed 
later. 



Figure 17 shows the contribution from each term of the — t- oscillation probability formula, 
Eq.(|3|). For sin^ 2^13 = 0.1(0.01) and 6 = 7r/2 with normal hierarchy, the contribution from the 



leading term, the CP violating (sin 5) term, and the matter term to the — )• z^e oscillation 
probability at 0.6 GeV neutrino energy are 0.05, -0.014, and 0.004 (0.005, -0.004, and 0.0004), 
respectively. The fraction of the contribution from the CP violating term to the total oscillation 
probability depends on the value of 6*13. For a relatively large value of ^13, the first term of Eq. ^ — 
which is proportional to sin^ ^13 — is dominant. Because the CP violating term is proportional to 
sin 6^13, for smaller 6*13 the asymmetry A = {P{i^^ z^e)-P(z7^ i7e)}/{P(f^ Ve)+P{T^ii T^e)] 
becomes larger. However, the number of — )■ Ve signal events decreases as Oiz gets smaller, 
resulting in larger statistical uncertainty. As a result, the sensitivity to the CP asymmetry has 
only modest dependence on ^13, if it is not too small to observe the signal. 

Due to the relatively short baseline and thus lower neutrino energy at the oscillation maximum, 
the contribution of the matter effect is smaller for the J-PARC to Hyper-Kamiokande experiment 
compared to other proposed experiments like LBNE in the United States [3]. Because the matter 
effect terms have sin^ 2^13 dependence, its significance relative to the sin (5 term (oc sin0i3) increases 
for a larger value of ^13. For sin^ 2^13 ~ 0.1, the matter effect has the same order of contribution 
to the oscillation probability as the CP violating term even with a baseline of 295 km, and thus 
the J-PARC to Hyper-Kamiokande experiment would have sensitivity to the mass hierarchy. 



2. J-PARC to Hyper-Kamiokande long baseline experiment 

The J-PARC neutrino beamline, which currently provides a neutrino beam to the T2K ex- 
periment, is designed to realize a common off- axis angle to Super-K and the candidate Hyper-K 
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FIG. 18. Expected neutrino flux at Hyper-K. Left: neutrino mode, right: anti-neutrino mode. i'^,i'^,i>(,,Ve 
components are shown separately. 



site in the Tochibora mine [M]. In this study, the off-axis angle is set to 2.5°, the same as the 
current T2K configuration. A beam power of 1.66 MW is assumed as the nominal case based on 
the KEK roadmap [55], while feasibility with lower beam intensity is also explored. The neutrino 
beam Monte Carlo simulation is used to estimate the neutrino flux and energy spectrum. In the 
beam simulation, the hadron production cross section is tuned based on the measurements of pion 
and kaon production cross sections by the NA61 collaboration |561 157) . A proton beam energy 



of 30 GeV and a magnetic horn current of 320 kA are assumed. Figure 18 shows the expected 
neutrino flux at Hyper-K for neutrino and anti-neutrino mode running. Thanks to the off-axis 
method, the spectrum has a narrow peak around the energy where oscillation probability is the 
maximum, with small high energy tail. Contamination of UeO^e) in the beam is well below 1% at 
the peak for both cases. 

Interactions of neutrinos in the Hyper-K detector is simulated with the NEUT program li- 
brary |58H60j ■ which is used in both the Super-K and T2K experiments. The response of the 
detector is simulated using the Super-K full Monte Carlo simulation based on the GEANT3 pack- 
age [61]. The simulation is based on the SK-IV configuration with the upgraded electronics and 
DAQ system, while the number of PMTs is reduced to about half to simulate 20% photocathode 
coverage of Hyper-K. Events are reconstructed with the Super-K reconstruction software. As de- 



scribed in Sec. II A , each Hyper-K tank is divided into five optically separated compartments, each 
one of which has about twice the volume of Super-K. Thus, the Super-K full simulation gives a 
realistic estimate of the Hyper-K performance. 

The electron (anti-)neutrinos from — )• fe oscillation are identified via their charged current 
interaction. The charged current quasielastic (CCQE) interactions, f g + n — )■ e~ + p and Ue+P ^ 



A Accelerator based neutrinos 



39 



e+ + n, have the largest cross section for E^, <\ GeV. Because the momentum of recoil protons are 
typically below Cherenkov threshold and neutrons do not emit Cherenkov light, they are identified 
as events with single electron-like ring in Hyper-K. Thanks to the two-body kinematics of the 
CCQE interactions, we can reconstruct the incident neutrino energy using only the energy and 
scattering angle of the outgoing electron(positron). Thus, we select events with a single electron- 
like ring as Vf, candidate signal events. The background is mainly intrinsic Ue contamination in the 
original neutrino beam, and neutral current interactions with 7r° misidentified as an electron due 
to mis-reconstruction of one 7 from tt*^ — >■ 77. 

The criteria to select f e candidate events are based on those developed for and established with 
the Super-K and T2K experiments: 

1. The event is fully contained (FC) inside the inner detector. 

2. Reconstructed vertex is inside the fiducial volume (FV). 

3. Visible energy (-Eyis) is greater than 100 MeV. 

4. Number of reconstructed rings is one. 

5. The reconstructed ring is identified as electron- like (e-like). 

6. There is no decay electron associated to the event. 

7. The invariant mass (Minv) of the reconstructed ring and a second ring force-found with a 
special 7r° fitter is less than 100 MeV/c^. This selection is imposed in order to reduce the 
background from mis-reconstructed -k^. 

8. The reconstructed energy {E^J'^) is less than 2 GeV. 

Assuming a charged current quasielastic interaction, the neutrino energy (El^^) is reconstructed 
from the electron energy (Eg) and the angle between the neutrino beam and the electron direction 
{6e) as 

^rec ^ 2(m„ - V)Ee + - {nin - Vy - ml 

2{mn-V -Ee+Pe cos Be) ' 

where mn,mp,me are the mass of neutron, proton, and electron, respectively, Pe is the electron 
momentum, and V is the nuclear potential energy (27 MeV). Here, rather than selecting a narrow 
signal region by cutting on the reconstructed energy, the reconstructed neutrino energy spectrum 
over a wider range is used to obtain maximum information on signal and background. Thus, the 
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TABLE VIII. The number of events after each selection step for neutrino mode, sin^ 26i3 =0.1, (5 = is 
assumed. 



1 / -1 r \ 

v mode (1.5 years) 


i^f, CC 


Vfi CC 


J^e CC 


CC 


NC 


CC 


J/^ l^e CC 


Interaction in FV 


64843 


2859 


4548 


368 


69898 


6024 


81 


FCFV & Evis>100 MeV 


46106 


2030 


3742 


297 


17437 


5858 


79 


1-ring 


24506 


1436 


2054 


180 


4201 


4972 


66 


e-hke 


711 


21 


1993 


175 


3211 


4893 


65 


No decay-e 


145 


4 


1639 


165 


2795 


4435 


64 


Mi„v <100 MeV/c2 


45 


1 


1173 


100 


811 


3966 


54 


< 2 GeV 


38 


1 


917 


57 


718 


3940 


51 


EfHciency(%) 


0.06 


0.03 


20.2 


15.4 


1.03 


65.4 


62.6 



selection criterion regarding reconstructed neutrino energy is much looser than the past studies. If 
we require 0.1 GeV< E^^^ < 1.25 GeV, the rejection efficiency for 'NClvr''' events, defined as neutral 
current events having no particle with momentum above the Cherenkov threshold other than a 
smg le vrO, is 95% while the efficiency for — z^g CC signal events is 64%. If necessary, a tighter 
selection can be applied by requiring Minv to be farther away from the vr^ mass and an additional 
requirement on the difference of likelihoods between the double and single ring assumptions {AC) 
as described in [46J. When Minv < 70 MeV/c^, AC < 20 and 0.1 GeV< E^J"" < 1.25 GeV are 
required, the NClvr'^ background rejection efficiency is 97.6% and the signal efficency is 50%. 



The number of events after each selection step is shown in Tables VIII and IX for each signal 
and background component. Unless otherwise stated, in this section running times of 1.5 years for 
neutrino mode and 3.5 years for anti- neutrino mode (five years in total) are assumed, with one year 
of running time corresponding to 10^ sec. For the signal, sin'^ 2^13 = 0.1 and 6 = are assumed in 
these tables. 



Figure 19 shows the reconstructed neutrino energy distribution after all the selections. In the 
neutrino mode, the dominant background components are intrinsic z^g contamination in the beam 
and mis-identified neutral current vr*^ production events originating from 1/^. In the anti- neutrino 
mode, in addition to I^g and 17^, and components have non- negligible contributions due to 
larger fluxes and cross-sections compared to their counterparts in the neutrino mode. 

The reconstructed neutrino energy distributions for several values of 6, with sin^ 2^13 = 0.1 and 



the normal mass hierarchy, are shown in the top plots of Fig. 20 The effect of 6 is clearly seen 
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TABLE IX. The number of events after each selection step for anti-ncutrmo mode, sin^ 2^13 ~ 0.1, (5 = is 
assumed. 



v mode [6.0 years) 


i^f, CC 


Vfi CC 


l^e CC 


Ve CC 


NC 


Vf. CC 


Ve CC 


Interaction in FV 


34446 


30121 


2905 


2335 


58569 


817 


2983 


FCFV & Evis>100 MeV 


24981 


20960 


2470 


1824 


17677 


798 


2907 


1-ring 


11830 


15925 


1168 


1225 


4433 


604 


2577 


e-hke 


401 


210 


1131 


1197 


3412 


596 


2539 


No decay-e 


88 


30 


893 


1148 


3000 


520 


2526 


Mi„v <100 MeV/c2 


20 


11 


564 


803 


834 


439 


2187 


< 2 GeV 


15 


10 


374 


598 


750 


421 


2168 


EfHciency(%) 


0.04 


0.03 


12.86 


25.63 


1.28 


51.6 


72.7 



600 



V mode 



400 F 



V mode 




FIG. 19. Reconstructed neutrino energy distribution of the candidate events, sin 2^13 = 0.1, J = and 
normal hierarchy. 



using the reconstructed neutrino energy. The bottom plots show the difference of reconstructed 
energy spectrum from (5 = for the cases 5 = jVTjTt and —\t^. The error bars correspond to the 
statistical uncertainty for the total five year experiment. By using not only the total number of 
events but also the reconstructed energy distribution, the sensitivity to 5 can be improved, and 
one can discriminate all the values of 5, including the difference between 5 = and tt. 

One has to note that, although the study described here is done within the standard MNS 
framework that attributes the source of CP violation to one parameter (5, the comparison of neu- 
trino and anti- neutrino oscillation probabilities gives a direct measurement of the CP asymmetry, 
without assuming an underlying model. 
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FIG. 20. Top: Reconstructed neutrino energy distribution for several values of S. sin^ 2013 = 0.1 and normal 
hierarchy is assumed. Bottom: Difference of the reconstructed neutrino energy distribution from the case 
with S — 0. The error bars represent the statistical uncertainties of each bin for 1.5 (3.5) years of running 
in neutrino (anti-neutrino) mode. 

3. Sensitivity to the CP asymmetry in the neutrino oscillation 

A analysis based on the reconstructed neutrino energy distribution has been performed to 
study the sensitivity of the ' J-PARC to Hyper-K' experiment to the CP asymmetry in the neutrino 
oscillation. 



Analysis method A binned is constructed from the E^^^ distribution, with 50 MeV bin 
width for the energy range of 0-2 GeV. As the systematic uncertainty, uncertainties in the normal- 
izations of signal, background originating from and i^^, those from and 17^, and the relative 
normalization between neutrino and anti-neutrino are taken into account. The is defined as 



X 



EE 



1 ± l^fu/l 



:i + /sig) • <g + (1 + /.J • <„ + (1 + 



Slg 



. 2 

Slg 



+ 



+ 



fl/e I fvjh 



where the index i runs over bins of reconstructed neutrino energy, and + and — are applied for 
neutrino and anti- neutrino mode, respectively. 
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A'^* is the number of expected events for the i-th E]^^ bin for oscillation parameters (^13"^, 5*'""*^, 
sign(Am32)*'^"'^). n^jg, nj,^, and n\,^ are the expected number of events for the appearance signal, 
the background originating from f^/v^, and the background from intrinsic v^jve^ respectively, with 
a set of oscillation parameter tested (^il'^*, 5*''^*, sign(Am32)*'^''''). The systematic parameters /sig, 
fvp.-, fvei fu/u represent uncertainties of the signal, the background from u^, those from i/g, and 
relative normalization of neutrino and anti-neutrino, respectively. Those systematic parameters 
are assumed to be energy independent and just the overall scale for each component is changed. 
ax is the assumed size of uncertainty for corresponding systematic parameter fx- 

In the recent search for appearance by T2K [1] , the total systematic uncertainty of the number 
of expected events at the far detector (Super-K) is 17.6% for sin^ 2^13 = 0.1, with contributions 
from neutrino flux uncertainty (8.5%), neutrino interaction cross section (10.5%), the near detector 
efficiency (5.6%), and the far detector related systematics (9.4%). 

The neutrino flux uncertainty is dominated by the hadron production uncertainty. The kaon 
production [57j and long target data [52] from the NA61 experiment will significantly reduce the flux 
uncertainty. Furthermore, systematic uncertainties from the neutrino beam and interaction cross 
section will be reduced by using energy dependent extrapolation from near detector measurements 
to far detector expectations. 

The neutrino interaction cross section and near detector systematics will be reduced with more 
data from the near detector. The T2K off-axis near detector consists of several subdetectors inside 
a 0.2 T magnetic field. In [Ij, inclusive i/^ charged current measurement in the near detector is 
used to normalize the expected event rate at the far detector, while the cross section uncertainty 
is estimated based on external measurements and interaction models. The uncertainty will be 
reduced by precise measurements of cross section with detailed analysis of the near detector data, 
which is in progress. 

For measurement of antineutrinos, background from the wrong sign component (i^ contamina- 
tion in P beam) introduces additional uncertainty if there is no charge sign selection capability. 
Thus, the magnetic field of the near detector will be a powerful tool to reduce uncertainty of 
anti-neutrino cross sections. Furthermore, the near detector may be upgraded for the J-PARC to 
Hyper-K experiment, based upon the experience with T2K. 

For the systematic uncertainties with the far detector, most of them are estimated by using 
atmospheric neutrinos as a control sample. For example, the uncertainty of the vr^ background 
rejection efficiency is estimated with the vr*^ topological control sample made by combining one data 
electron and one simulated gamma event; its uncertainty is limited by the number of atmospheric 



44 III PHYSICS POTENTIAL 

Hyper-K (560kt FV) / 1 .5yrs v + 3.5yrs v /1 .66MW 



S 

to 



'0 0.02 0.04 0.06 0.08 0.1 0.12 
sin22ei3 

FIG. 21. Allowed regions for the parameter sets of sin^ 2^13 = (0.01,0.03,0.05,0.07,0.10) and S = 
(— ^TT, 0, ^TT, tt) overlaid together. Blue, green, and red lines represent 1, 2, 3 tr allowed regions, respec- 
tively. Stars indicate the true parameters. It is assumed that the mass hierarchy is known to be the normal 
hierarchy. 

neutrino events available. Similarly, ring counting and particle identification uncertainties are lim- 
ited by statistics of the z^e-enriched atmospheric neutrino sample. Another source of uncertainty 
is the limited knowledge of the neutrino interaction cross section. With more than an order of 
magnitude larger statistics available with Hyper-K, the beam neutrino events as well as the atmo- 
spheric neutrino events can be used to study systematics. Together with improved understanding 
of the neutrino interaction, uncertainties associated with the far detector will be reduced. 

Based on the experience from T2K analysis and prospects for future improvements described 
above, ax is set to 5% for all four systematic parameters. 

For each set of (^il*^*, 6^°^^, sign(Am32)*^'^*), the is minimized by changing the systematic 
parameters, fx- The is then compared to the value at the true oscillation parameters, and the 
difference Ax^ = (x^(test) — x^(true)) is used to evaluate the significance of the measurement. 
When allowed regions are drawn on a (sin^ 26*13 )-5 plane, the 68.3% (la), 95.5%(2a), and 99.7% 
(3(t) CL allowed regions are defined as the regions of parameters where Ax^ < 2.30,6.18, and 
11.83, respectively. For measurements where a single parameter is concerned, e.g. the uncertainty 
of 6, mass hierarchy determination, and exclusion of sinS = 0, Ax^ values corresponding to a single 
parameter are used [5]. 

Sensitivity if the mass hierarchy is known Let us first consider the case where the mass 
hierarchy is known by other measurements. If the mass hierarchy is known to be the normal 
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FIG. 22. Allowed regions for inverted hierarchy. See caption of Fig. [2T| 
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FIG. 23. 1(7 error of (5 as a function of sin 26'i3 for the normal hierarchy case. 
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FIG. 24. Sensitivity to GP violation. Blue, green, and red lines correspond to 1, 2, and 3 a exclusion of 
sin (5 — 0, respectively. Left: horizontal axis is linear scale, right: log scale. 
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FIG. 25. Reconstructed neutrino energy distributions for (sin^ 2^13, (5) = (0.10, O.Stt) with normal mass 
hierarchy (black solid) and (sin^ 26*13, (5) = (0.10, 0.9n) with inverted mass hierarchy (red dashed). 

hierarchy, sign(Am32) can be fixed to be positive. The allowed regions for the true parameter sets 
of combinations of sin 26*13 = (0.01,0.03,0.05,0.07,0.1) and 6 = (0, ^vr, vr, -^vr), for the normal 
hierarchy, are shown in Fig. [21] The same plot for the inverted hierarchy case is shown in Fig. [22} 
For both cases, 6 and sin'^ 2^13 can be well determined for the region of sin^ 2^13 indicated by the 
T2K result of June 2011 [Ij, sin^ 26*13 > 0.03(0.04) (90% CL) for normal (inverted) hierarchy. 
Figure 23 shows the la error of (5 as a function of sin^ 2^13 for the normal hierarchy. The size 



of the error is almost the same for the inverted hierarchy. For sin^ 2^13 > 0.03, the value of 6 can 
be determined to better than 8° for 5 = 0°, and better than 18° for 5 = 90°. 



Figure 24 shows the regions where sin 5 = is excluded, i.e. CP is found to be violated in 



the lepton sector, with la, 2a, and 3a significance, in the case of normal hierarchy and the mass 
hierarchy is known. The J-PARC Hyper-K experiment will have a sensitivity to the CP asymmetry 
down to sin'^ 2^13 ~ 0.003 with 3 a significance. 



Effect of unknown mass hierarchy If the mass hierarchy is unknown, there might be a fake 



solution due to the parameter degeneracy [63]. Figure 25 shows the reconstructed neutrino energy 
distributions for (sin^ 2^13, 6) = (0.10, 0.57r) with normal mass hierarchy, and (sin^ 2^13, 6) = (0.10, 
0.97r) with inverted mass hierarchy. Those two parameters give similar oscillation probabilities, 
hence similar energy distributions. 

Figure 26 shows Ax^ as a function of 5, for the case the mass hierarchy is unknown. In this plot. 



sin^ 2^13 is fixed to 0.1. The true parameters are (sin^ 2^13, 6) = (0.10, 0.57r) and normal hierarchy. 
In the case the mass hierarchy is unknown, is calculated with both normal and inverted hierarchy 
hypotheses and the smaller is taken to calculate A^^- Due to the parameter degeneracy, there is 
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5[Ji] 

FIG. 26. The Ax^ as a function of S, for the case the mass hierarchy is unknown, sin^ 2^13 is fixed to 0.1. 
The true parameters are (sin^ 2^13, 5) = (0.10, O.Stt) and normal hierarchy. In this case, is calculated with 
both normal and inverted hierarchy and the smaller is taken as the minimum value. Due to the parameter 
degeneracy, there is a local minimum with inverted hierarchy hypothesis. 



a local minimum with the inverted hierarchy hypothesis in addition to the minimum around the 
true value with the normal hierarchy hypothesis. If this fake solution is consistent with null CP 
asymmetry even though the true solution violates the CP symmetry, then the sensitivity to the 
CP violation will be lost for that parameter set. 



Figure 27 shows the 3a allowed regions for the case the true mass hierarchy is normal but 
not determined prior to this experiment. Solid red (dashed blue) line shows the contour for true 
6 = 0(271"). Stars indicate the true parameter values. Plots for sin^ 2^13 of 0.01, 0.05, and 0.1 are 
shown together. The areas indicated by dashed lines around sin^ 2^13 = 0.1 correspond to the case 



shown in Fig. 26 In addition to the region around the true values, there is a fake solution due to 
unknown mass hierarchy. 

If sin^ 2^13 is as large as ~ 0.1, the NOi^A experiment in United States [17] together with 
the T2K experiment will have sensitivity to the mass hierarchy. In addition, for a large value of 



sin^20i3, Hyper-K itself has considerable sensitivity to the mass hierarchy. Figure 28 shows the 
sensitivity to the mass hierarchy from the J-PARC to Hyper-K experiment alone. For each set of 
(0*3^^, S^'^^^) with normal mass hierarchy (Am^2^'^^^ > 0), the inverted hierarchy hypothesis is tested 
by calculating Ax^ for various sets of {e\f\ Jt^*^*) with Aml^ < 0- The J-PARC to Hyper-K 
experiment has sensitivity to the mass hierarchy for sin^ 2^13 > 0.05. For sin^ 2^13 = 0.1, the mass 
hierarchy can be determined with more than 3a significance for 46% of the 6 parameter space. Also, 
the atmospheric neutrino observation will have a sensitivity to the mass hierarchy as described in 
Sec. IIIB If sin^ 2^13 is too small for those experiments to determine the mass hierarchy, the effect 
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FIG. 27. 3(T allowed regions for the case the mass hierarchy is unknown. The true mass hierarchy is 
the normal hierarchy. Solid (dashed) line shows the contour for true S — 0{-^n). Stars indicate the true 
parameter values. Plots with sin^ 26*13 of 0.01, 0.05, and 0.1 are overlaid. In addition to the region around 
the true values, there are fake solutions around S = 7r(— tt) for true 6 = 0(^7r) due to degeneracy. 



Hyper-K (560kt FV) / 1 .5yrs v + 3.5yrs v/1 .66IVIW 





0.05 0.1 
Sin22ei3 



FIG. 28. Sensitivity to the mass hierarchy. Blue, green, and red lines correspond to 1, 2, and 3 a, respectively. 
The true mass hierarchy is the normal hierarchy. 



of mass hierarchy becomes smaller, as it becomes relatively insignificant compared to the genuine 
CP asymmetry. 

If we assume that the mass hierarchy is not known, the discovery potential for CP violation 



is as shown in Figure 29 For the parameter < 5 < vr, the exclusion region becomes smaller 



compared to the case where mass hierarchy is known (Fig. 24 ) because of the fake solution coming 



from unknown mass hierarchy. However, there are sets of parameters for which CP asymmetry can 
be observed, and for — tt < 5 < 0, there is little effect even if the mass hierarchy is not measured. 




FIG. 29. Sensitivity to CP violation for the case the mass hierarchy is not known. Blue, green, and red 
lines correspond to 1, 2, and 3 a exclusion of sin 5 — 0, respectively. The true mass hierarchy is the normal 
hierarchy. 




FIG. 30. Fraction of b for which sin (5 = is excluded with 3 a significance as a function of true sin^ 20i3, 
with different ratio of v and v running time, while total running time is fixed to five years with 1.66 MW. 
Normal hierarchy is assumed. 



This is because the mass hierarchy can be determined and the degeneracy is resolved for these 



parameter sets as shown in Fig. 28 For the case in which the true mass hierarchy is inverted, a 
similar argument holds with — vr < 5 < and < (5 < vr inverted. 

In reality, the mass hierarchy can be determined for wider range of parameters by combining 
other experiments and atmospheric neutrino observations with Hyper-K. In conclusion, the knowl- 
edge of the neutrino mass hierarchy will have only limited impact on the discovery potential of 
leptonic CP violation by a J-PARC to Hyper-K long baseline program as discussed here. 
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FIG. 31. 1(7 uncertainty of (5 as a function of integrated beam power. The ratio of v and v running time is 
fixed to 7:3. left: sin^ 26*13 = 0.1, right: sin^ 26*13 = 0.01. 
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FIG. 32. Fraction of 5 for which sin (5 — can be excluded with 3 cr as a function of the integrated beam 
power. The ratio of neutrino and anti-neutrino mode is fixed to 3:7. 



Sensitivity vs running time Figure 30 shows the fraction of 5 for which sin 5 = is excluded 
with 3 a, with different ratio of v and v' running time, while the total running time and the beam 
power is fixed to five years with 1.66 MW. The case with v : = 3 : 7 has the best sensitivity, 
although the difference is marginal around this value. 



Figure 31 shows the la uncertainty of 5 as a function of the integrated beam power for sin^ 2^13 = 
0.1 and 0.01. The mass hierarchy is assumed to be known. The ratio of v and v' running time is 
fixed to 7:3. 



Figure 32 shows the sensitivity to CP violation as a function of the integrated beam power. The 
vertical axis shows the fraction of 6 for which sin 5 = is excluded with 3a significance. The ratio 
of neutrino and anti-neutrino mode is fixed to 3:7. Solid and dashed lines correspond to the case 
the mass hierarchy is known and unknown, respectively. The true mass hierarchy is normal in both 
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FIG. 33. Fraction of d for which sinS = is excluded with three sigma as a function of the true value of 
sin^ 2^13. Red and Black lines represent 3 MW-year and 8.3 MW-year, respectively. 



cases. Although the sensitivity becomes worse if the mass hierarchy is unknown, CP asymmetry 
can be accessed for more than half of 6 parameter region if sin^ 2^13 > 0.01. 



Figure 33 shows the fraction of S for which sin 5 = is excluded with 3a as a function of the 
true value of sin^20i3, for 3 MWxlO'^s and 8.3 MWxlO^s of integrated beam power. The mass 
hierarchy is assumed to be normal hierarchy and measured prior to this experiment. 

The fraction of 6 for which CP asymmetry can be discovered with more than 3a under several 
assumptions are summarized in Table IXl 



TABLE X. Fraction of S (in %) for which CP asymmetry can be discovered with more than 3a under several 
assumptions. Integrated beam power is 8.3 MW-year. 





mass hierarchy 




known 


unknown 


0.1 


74 


55 


sin^ 26*13 0.03 


74 


63 


0.01 


66 


59 



4-. Summary 

The sensitivity to leptonic CP asymmetry of a long baseline experiment using a neutrino beam 
directed from J-PARC to the Hyper-Kamiokande detector has been studied. The running time is 
assumed to be five years in total, with 1.66 MW of beam power. Based on a full simulation of 
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beamline and detector, it is found that S and sin^ 2^13 can be well determined for the values of 
sin^ 2^13 indicated by the recent T2K result. 

If the mass hierarchy is known, for sin^ 2^13 > 0.03 the value of 6 can be determined to better 
than 18° for all values of S and CP violation in the lepton sector can be observed with 3a significance 
for 74% of the possible values of 5. If we assume that the mass hierarchy is not known, the sensitivity 
to CP violation is reduced due to degeneracy. Even for this case, CP violation can be observed 
with 3(7 significance for 55% of S parameter space if sin^ 2^13 = 0.1, and 63% if sin^ 2^13 = 0.03 
with five years of experimental operations. 

For sin^ 2^13 > 0.05, it is also possible to determine the mass hierarchy for some of S with the 
J- PARC to Hyper-K experiment alone. For sin^ 2^13 = 0.1, the mass hierarchy can be determined 
with more than 3a significance for 46% of the S parameter space. 
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B. Atmospheric neutrinos 

1. Goals of the atmospheric neutrino study 

Atmospheric neutrinos are a guaranteed neutrino source in the Hyper-Kamiokande experiment. 
The indication by T2K [Ij that ^13 is potentially large is encouraging for future atmospheric 
neutrino studies, as there will be a good chance to extract information on neutrino properties via 
the three flavor oscillation effect. Assuming sin^ 2^13 > 0.04 as the global fit result suggests [28], 
the targets of the atmospheric neutrino studies in Hyper-K would be: 

• mass hierarchy determination, namely to select > or < with more than 3a 
significance provided sin^ 623 > 0.4. 

• to solve sin^ 623 octant degeneracy, namely to discriminate sin^ 623 < 0.5 (first octant) from 
sin^ 923 > 0.5 (second octant), when the mixing is not maximal as sin^ 2^23 < 0.99. 

• to obtain complementary information on CP phase 6. 

To extract the expected three flavor oscillation effects, we will study atmospheric electron neutrino 
flux variations as well as muon neutrino flux variations. Expected sensitivities for all of these topics 
are discussed in this section. 

Oscillation probabilities of atmospheric neutrinos in the three flavor neutrino mixing scheme 
have been discussed by many authors, and the oscillation effect in electron neutrino flux is analyt- 
ically calculated [64J as: 



where we call the first, second, and third terms the "solar term", "interference term", and "^13 
resonance term", respectively. P2 is the two neutrino transition probability of i^e — t- f^^,- which is 
driven by the solar neutrino mass difference Am2i. R2 and I2 represent oscillation amplitudes for 
CP even and odd terms. For anti- neutrinos, the probabilities P2,R2,l2 are obtained by replacing 
the matter potential V — )• —V, and the sign of the 6 (see f64j for details), r is the v^/fe flux ratio 
as a function of neutrino energy; r ~ 2 at sub-GeV energies, starts deviating from 2 at 1 GeV, and 
reaches to ~ 3 at 10 GeV. The ^13 is an effective mixing angle in the Earth; sin^ ^13 could become 
large at 5 ~ 10 GeV neutrino energy due to the matter potential [65ti67] . This MSW resonance 



1 « P2 • (r • COS^ 6*23 - 1) 

— r • sin ^13 • cos^ 613 ■ sin 26*23 ■ (cos 5 ■ R2 — sin 5 ■ I2) 
+2sin2 ^i3-(r-sin2 ^23-!) 



(5) 
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FIG. 34. Oscillated i^e flux relative to the non-oscillated flux as a function of neutrino energy for the 
upward-going neutrinos with zenith angle cos0,y = —0.8. Pe is not included in the plots. Thin solid 
lines, dashed lines, and dotted lines correspond to the solar term, the interference term, and the ^13 
resonance term, respectively (see Eq. [s]). Thick solid lines are total fluxes. Parameters are set as 
(sin^ 012, sin^ 6I13, sin^ 6*23, S, Am^i, Am§2) = (0.31, 0.025, 0.6, 40°, 7.6 x 10~^eY^, -f 2.4 x lO-^eV^) unless oth- 
erwise noted. The 623 octant effect can be seen by comparing (a) (sin^ 6*23 = 0.4) and (b) (sin^ 6*23 = 0.6). S 
value is changed to 220° in (c) to be compared with 40° in (b). The mass hierarchy is inverted only in (d) 
so 6*13 resonance (MSW) effect disappears in this plot. For the inverted hierarchy the MSW effect should 
appear in the Pe flux, which is not shown in the plot. 



happens with neutrinos in the case of normal mass hierarchy > 0)) ^^^d with anti- neutrinos 

in the case of inverted mass hierarchy < 0)- 



In order to demonstrate the behavior of these three terms. Fig. 34 shows how the fe flux changes 
as a function of neutrino energy based on a numerical calculation of oscillation probabilities, in 
which the matter density profile in the Earth is taken into account \25\ I68j . We adopted an 
Earth model constructed by the median density in each of the dominant regions of the preliminary 
reference Earth model (PREM) [69]: inner core (0 < r < 1220km) 13.0 g/cm^, outer core (1220 < 
r < 3480km) 11.3 g/cm^, mantle (3480 < r < 5701km) 5.0 g/cm^, and the crust (5701 < r < 



6371km) 3.3 g/cm^. In Fig. 34 dotted lines correspond to the ^13 resonance term (the third term 
in Eq. 5^, which could make a significant contribution in the 5 ~ 10 GeV region if sin^ ^13 is a few 
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percent [I1[2H]. The resonance effect in (not Ve) can be seen in Fig. 34 ^b) which assumes normal 



mass hierarchy, but resonance does not occur in Fig. 34 d) which assumes inverted hierarchy. 
For the case of i^e, the situation is reversed. By using statistically enhanced Ve and samples 
in observed data, this resonance term would enable us to reveal the neutrino mass hierarchy. In 
addition, this term (as well as the solar term) has discrimination power for the 623 octant because 



this term is proportional to (r • sin^ 623 — 1). The effect of 623 can be seen by comparing Fig. 34 a) 
and (b). Dashed lines in Fig. |34] correspond to the interference term (second term in Eq. [5]). The 
CP phase effect would be dominant in the region from a few 100 MeV to a few GeV neutrino 
energy as shown in Fig. [34|b) [6 = 40°) and (c) {6 = 220°). This interference term is proportional 
to sin 2^23 (the first power) but the ambiguity of the parameter is as small as a few % and therefore 
CP sensitivity does not much depend on the value of 023- 



2. Analysis method 

Atmospheric neutrino interactions are simulated by the Monte Carlo method using a flux calcu- 
lation [To] and interaction models |58H60j . The propagation of secondary particles and Cherenkov 
photons in detector water and event reconstructions are performed by using the Super-K detector 
simulator and standard reconstruction tools [20] . We expect that - by design - the detector per- 
formance of Hyper-K will be basically the same as that of Super-K. In addition to using similar 
photodetectors and electronics, and having similar quality water, the length scale of each seg- 
mented Hyper-K subdetector compartment is 50 m, which is only 1.5 times larger than that of 
Super-K. We have generated Monte Carlo (MC) events with 11 Megaton-years statistics in the 
fiducial volume which correspond to a 500 year exposure of Super-K or 25 years of Hyper-K. Os- 
cillation probabilities are calculated by taking into account the full parameters in the standard 



3 flavor neutrino scheme; ^12, ^13, 023, and two squared mass differences of Am^i and Am? 



32- 



Oscillation parameters, if they are not otherwise specifled, are set to the values shown in Table XI 
The ^23 value is set to be the maximal mixing suggested by atmospheric and accelerator neutrino 
experiments |20H23t [52] , and other parameters are taken from a recent global fit result [28] except 
for the random choice of the 6 value. In the mass hierarchy determination study, we test both signs 
of Am'^2- 

The MC events are divided into three classes: fully contained (FC), partially contained (PC), 
and upward- going muons (UP;u). The FC events are further divided into several sub-samples 
based on the Super-K analyses by using reconstructed variables such as visible energy, the number 
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TABLE XI. Oscillation parameters are set to be these values in this section, unless otherwise noted. 



A 2 


7.6x10 ^ eV 


ATO32 


2.4x10-3 eV^ 


sin^ 6*12 


0.31 


sin^ 013 


0.025 


sin^ 6*23 


0.5 


6 


40° 



TABLE XIL Expected number of I'e-like and i^e-\ike events in 10 Hyper-K years for each interaction com- 
ponent. 





CC l^e 


CC De 


CC + 


NC 


Total 


i^e-like sample 


15247 


2831 


3731 


4792 


26601 


- single-ring 


6356 


1086 


1682 


1740 


10864 


- multi-ring 


8891 


1745 


2049 


3052 


15737 


Percentage (%) 


57.3 


10.6 


14.0 


18.0 


100.0 


Pg-like sample 


28309 


17255 


1232 


4559 


51355 


- single-ring 


20470 


13401 


444 


2496 


36811 


- multi-ring 


7839 


3854 


788 


2063 


14544 


Percentage (%) 


55.1 


33.6 


2.4 


8.9 


100.0 



of Cherenkov rings, particle type (e-like or //-like), the number of muon decay electrons, and so on 
[25]. In addition, the multi-GeV single-ring and multi-ring e-like events are regrouped into fg-like 
and Pg-like samples in the following manner. In the charged current (CC) non quasi-elastic (QE) 
interaction components among the single-ring sample, 7r+(7r-) is expected to be more copiously 
produced in t'e(^'e) interactions because the secondary charged lepton is e~{e'^). The vr"'" decays 
into n'^ which in turn produces a delayed signal of a decay electron, while in the case of tt~ 
it is often absorbed in water before decaying into /x^, so no decay electron would be produced. 
Therefore, events with more than one decay electron are classified as fg-like, while events having 
no decay electron are classified as z^e-like. As for the multi-GeV multi-ring e-like sample, the events 
are also divided into fg-like and z^g-like samples via a likelihood method. CC Ve interactions tend 
to have a larger Feynman y distribution than CC Ug, therefore CC z^e is expected to have larger 
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transverse momentum, more rings, and more muon decay electrons. Hence these three observables 
are used in the construction of a hkehhood function, and multi-ring events are divided into i/g-like 



and z^e-like sub-samples by applying the likelihood cut. Table XII shows the expected number of 
non-oscillated fe and Pe-hke events in 10 years exposure (for the z^e-like sample: CC z^e 57.3%, CC 
z^e 10.6%, while for the i^g-like sample: CC Ve 55.1%, CC i>e 33.6%). Due to the small fraction of 
wrong sign electrons (CC z^e ~ 10%) in the fe-like sample, we expect a significant enhancement 
of the CC Ve component (CC fe ~ 60%) in the case of the normal hierarchy, while less excess is 
to be expected in the case of the inverted hierarchy. The i/g-like sample will provide independent 
information on the mass hierarchy with its higher CC z^e fraction. 

These MC events are used as fake data as well as expectation. To make fake data, the MC is 
weighted by livetime and oscillation probabilities. The is defined by comparing each fake data 
set with the expectation in terms of Poisson statistics: 



+ E(^) (6) 



where 



n : counter for event type, momentum, and zenith angle 
: the number of MC events in the n-th bin 
A^p-p : the number of fake data events in the n-th bin 
i : counter for uncertainties 

€i : the list of systematic uncertainties (nuisance parameters) 

: error coefficients for i-th uncertainty and n-th event bin 
(Tj : estimated size of systematic uncertainties 



The definition of and data binning are same as the latest Super-K three fiavor neutrino oscillation 
analyses [25]. Moreover, systematic uncertainties (fij) based on [25] which covers uncertainties of 
neutrino fiux, interactions, detector response, and event reconstructions are incorporated in this 
analysis. Because the fake data are made from MC events, minimum is always obtained as 
Xmin ~ with all ej = at the true oscillation parameters. Then sensitivities of each oscillation 
parameters are evaluated by using Ax^ = Xmin(t6st point) — Xmm(true point) = Xmin(test point). 
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FIG. 35. Expected event rate changes in (a) sub-GeV single-ring e-like, (b) multi-GeV z^e-like, and (c) 
niulti-GeV i?e-like event samples. The vertical axis shows the ratio of oscillated e-like event rate to the 
non-oscillated one. Mass hierarchy is normal for dashed lines and inverted for solid lines. Colors show 
sin^ 023 values as 0.4 (black), 0.5 (green), and 0.6 (red). Points with error bars represent null oscillation 
expectations with expected statistical errors for 5.6 Megaton-years exposure or 10 years of Hyper-K. 



3. Neutrino mass hierarchy 



In Fig. 35, the expected zenith angle distributions of e-like events are shown separately for 
sub-GeV e-like, multi-GeV fg-like, and multi-GeV z/e-like sub-samples as the ratio against the non- 
oscillated case. We expect a sizable difference between normal (dashed lines) and inverted (solid 
lines) hierarchy both in the z^e-like and z^g-like samples. The difference becomes larger for larger 
sin^ 023 because the resonance term is proportional to (r • sin'^ 623 — 1) in Eq. Isl Sensitivity for 



determining the neutrino mass hierarchy is studied and shown in Fig. 36 The Ax for the wrong 
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Hyper-K years Hyper-K years 



FIG. 36. Expected significance for mass hierarchy determination as a function of Hyper-K hvetime in years. 
All of ^23: ^13: and 6 are assumed to be unknown and allowed to vary without any constraints. 




Hyper-K years Hyper-K years 



FIG. 37. Expected significance for mass hierarchy determination as a function of Hyper-K livetime in years. 
Both 6*23 and ^13 are assumed to be well known and are fixed. 

mass hierarchy assumption, - interpreted as significance of the mass hierarchy determination - is 
shown as a fmiction of livetime in years. In the calculation of the significance, fitting parameters 
^13) ^23) and 6 are assumed to be unknown and allowed to freely vary when obtaining the 
minimum. The expected significance is larger for larger sin^ 623 and sin^ ^13 as is expected. With 
10 years exposure, the significance is more than 3a for most of the parameter sets suggested by the 



global fit results [28] . Figure 37 also shows the significance but in this case both ^13 and 6*23 are 
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FIG. 38. Expected significance for the mass hierarchy determination. In the left panel, normal mass 
hierarchy is the case and for the wrong assumption; Ax^ = x^^^^{mverted) — x^j^^ (normal) is shown for 
various true values of sin^ 2^13. The right panel is for the inverted hierarchy case. Each colors show the case 
of sin^6'23 — 0.4 (black), 0.5 (red), and 0.6 (green), and the blue horizontal lines show 3(t(Ax^ ~ 9.2). 

assumed to be well known and are fixed. The effect of the assumption that both mixing angles are 
known or not can be seen most readily in the inverted hierarchy case, but even there the difference 
is rather hmited. In the case of (sin^ 26*13, sin^ 6*23) = (0.08,0.5) for example, required exposure 
time to confirm inverted hierarchy with 3cr CL would be between 5 years (^13 and 623 are well 
known) and 8 years (^13 and 623 are unknown). On the other hand, we need ~ 5 years to confirm 
normal hierarchy no matter what the state of knowledge regarding ^13 and 623 is. In conclusion 
we need between 2 ~ 10 years of data to reach 3a significance for the parameter space suggested 
by the global fit results except for the case that sin^ 623 is smaller than 0.4. Figure 
the 10 year significance for each sin^ 623 and sin^ 613 parameter set. 
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summarizes 



4-. Octant of 9 



23 



Figure 35 already showed how the event rate changes for different sin^ 623 values. We expect 
event rate change and zenith angle distortions in multi-GeV fe and z/e-like samples as well as small 
rate changes in sub-GeV electrons. 



Figure 39 shows the expected sensitivities for sin^ 623 and sin^ ^13 with an exposure of 10 Hyper- 
K years. Normal mass hierarchy is assumed in the figure. We can expect to discriminate between 
sin^ 623 = 0.4 and 0.6, which corresponds to sin^ 26*23 = 0.96. In the case of sin^ 623 = 0.45 or 0.55, 
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FIG. 39. Expected sensitivities for sin^ 6*23 and sin^ 26*13 at 90% CL with a livetime of 10 years in Hyper- 
K. Stars in the contours represent the assumed true mixing angles. The left panel shows the case of 
sin^ 26*23 = 0.96 (sin^ 6*23 = 0.4 or 0.6). The right panel shows the case of sin^ 26*33 = 0.99 (sin^ 6*23 = 0.45 
or 0.55). Normal mass hierarchy is assumed. 

which corresponds to sin^ 2^23 = 0.99, discrimination is marginal but could be achieved. 



5. CP phase 

The effect of CP phase 5 expressed by the second term in Eq. [5] can appear in the neutrino 



energy region from 100 MeV to a few GeV as shown in Fig. [34} Figure 40 shows that the event 
rate changes will appear both in sub-GeV and multi-GeV energy samples as expected. 

Figure 41 shows the expected sensitivity for the CP phase evaluated by Ax^. Normal mass 



hierarchy is assumed in this analysis. These sensitivities do not depend on 623 very much and 
so sin^ ^23 is fixed to be 0.5 in the figure. In the case of sin^ 2^13 = 0.16, discrimination of the 
CP phase for 5 = 40°, 140°, 220°, and 320° is expected to be possible at the 90% CL and even 
99% CL discrimination is possible for S = 40° and 220°. In the case of sin^ 2^13 = 0.08 the 
accuracy of CP phase measurements gets worse, but constraints could still be obtained in some 
cases. Although the sensitive parameter space and 5 precision are limited, there is a good chance 
to obtain precious information on CP 6 and to provide supplemental information for the CP 
asymmetry study conducted by the J-PARC to Hyper-K neutrino oscillation experiment discussed 
in Sec. HITAl 
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FIG. 40. Expected event rate changes in sub-GeV single-ring e-like, multi-GeV t/g-like, and multi-GeV 
Pg-like event samples. The vertical axis shows the ratio of oscillated e-like event rate to the non-oscillated 
one. True CP phase S is varied as 40° (black), 140° (red), 220° (green), and 320° (blue). Normal mass 
hierarchy is assumed. Points with error bars represent null oscillation expectation with expected statistical 
errors for 10 years of Hyper-K. 

6. Summary 



Sensitivity to several neutrino oscillation parameters by using high statistic atmospheric neu- 
trino data in Hyper-Kamiokande has been studied assuming sin^ 2^13 > 0.04. We have applied 
the full MC simulation and latest systematic uncertainties used in the current Super-Kamiokande 
analyses. It has been found that the expected significance for the mass hierarchy determination is 
more than 3a provided sin^ 623 > 0.4. We expect to be able to discriminate between sin^ 623 < 0.5 
(first octant) and > 0.5 (second octant) if sin^ 2^23 is less than 0.99. In addition, there is a good 
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FIG. 41. Expected sensitivities for 6 and sin^ 26*13 at 90% CL (red) and 99% CL (blue) with a livetime of 10 
Hyper-K years. Stars in the contours represent the true parameter set of ^13 and 6. Normal mass hierarchy 
is assumed. 



chance to obtain supplemental information on CP 5. 
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C. Nucleon decays 

Large water Cherenkov detectors have very good sensitivities for nucleon decays; indeed, the 
technology was originally developed and the first generation of such detectors was constructed in the 
1980's primarily to search for proton decays predicted by the then-popular Grand Unified Theory 
(GUT), Minimal SU{5). That model was eventually ruled out through non-observation of decays, 
leading to new models with longer lifetime predictions based on the experimentally demonstrated 
rarity of these decays. The search for these predicted decays continues: for more than a decade 
Super-Kamiokande has had the world's best limits, generally by an order of magnitude or more. 



on most of the current theoretically favored decay modes. As described in Sec. IA2 among many 
possible nucleon decay modes p — t- e'^n^ and p — t- uK'^ have been the subjects of the most intense 
interest, and they will be discussed in this section. 

The sensitivity of the Hyper-Kamiokande experiment for nucleon decays has been studied with a 
MC simulation based on the Super-Kamiokande analysis. An estimate of the atmospheric neutrino 
background is necessarily included in the study. 



1. Sensitivity study for the p — ?> c+tt'^ mode 

For the p — )• e'^n^ mode where a vr^ decays primarily to two 7s, all of the final state particles are 
visible in a water Cherenkov detector. Signal candidates are selected with the following criteria: 
(A-1) the event is fully contained in the detector (FC event), and the vertex is reconstructed in 
the fiducial volume of 0.56 Megatons, (A-2) all rings are e-like, and the number of rings is two 
or three (two ring events could occur due to an overlap of rings), (A-3) there is no Michel decay 
electron (eliminate the events with an invisible muon) , ( A-4) for three ring events only, the invariant 
mass of two e-like rings is reconstructed between 85 to 185 MeV/c^ to identify the vr'^, (A-5) the 
reconstructed total momentum is less than 250 MeV/c and the total invariant mass is between 800 
to 1050 MeV/c2. 

Water molecules consist of two free protons and eight bound protons in an oxygen nucleus. In 
the MC simulation, the binding energy, Fermi motion, and interactions in nucleus are taken into 
account. 

A significant but unavoidable inefficiency of the p — )• e~^TT^ signal is caused by missing the 
which interacts in the nucleus (absorption and scattering). Although the p — )• e^n^ efficiency of 
a free proton is as high as 87%, the average efficiency of both free and bound protons is lower 
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because the vr*^ from the bound proton decay often suffers these interactions inside the nucleus. As 
a result, the overall proton decay efficiency of p — t- e+vr'' is estimated to be 45%. 

The main background source of the proton decay search is the atmospheric neutrino events, 
which can occasionally produce an electron and a vr^ in the final state. With the selection criteria 
of (A-l)-(A-5), the remaining background events are estimated to be 1.6 events/Megaton-year from 
the atmospheric neutrino MC simulation. This result of the MC simulation has been experimentally 
confirmed by the K2K experiment [7l]. K2K's one kiloton water Cherenkov near detector accu- 
mulated muon neutrino beam data corresponding to a 10 Megaton-year exposure of atmospheric 
neutrinos. Using a calibrated and well-understood neutrino fiux, the corresponding atmospheric 
neutrino-induced background rate of false p — )• fi^ir^ events was estimated based on the observed 
beam-induced events in the detector. This rate was then extrapolated in turn to the atmospheric 
neutrino-induced background rate of p — >■ e'^n^, 1 .631 o;33 (stat)t (syst) events /Megaton •year. 

Fig. [42] shows the total invariant mass distributions of the proton decay signal events and the 
atmospheric neutrino background events after all the cuts except the mass cut. The running time of 
Hyper-Kamiokande is assumed to be 10 years which corresponds to a 5.6 Megaton- year exposure. 
The proton lifetimes in these plots are taken to be: (a) 1.2x10'^^ years (current limit), (b) 2.5x10^^ 
years, (c) 5.0x10^'^ years, and (d) 1.0x10^^ years. The number of atmospheric neutrino background 



events is estimated to be 9.0 events. Fig. 43 shows the significance of the signal events as a function 
of the proton lifetime by counting the events in the signal region. In this figure, the systematic 
uncertainty of the background is assumed to be 44%, the same as Super-K. We could separate a 
signal from the background events at a 5 o" significance for a proton lifetime of 3.4x10^^ years, and 
at 3 fj with 5.7x10^^ years. 



Fig. 44 shows the sensitivity for proton decay with a 90% CL as a function of the detector 
exposure. We could reach 1.0x10^^ years partial lifetime with a 4 Megaton-year exposure, which 
corresponds to 8 years running Hyper-Kamiokande; by contrast, it will take Super-Kamiokande 
178 years to reach this level. 



2. Sensitivity study for the p — >■ h'K'^ mode 

For the p — t- uK'^ mode, K'^ itself is not visible in a water Cherenkov detector due to having a 
low, sub-Cherenkov threshold, momentum. However, can be identified by the decay products 
of —7- ^+ + v (64% branching fraction) and — )• 7r+ + (21% branching fraction). The 
muons and pions from the decays have monochromatic momenta due to being produced via 
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FIG. 42. Invariant mass distributions of MC events after all cuts except for the mass cut with a 
5.6 Megaton-year exposure. The proton lifetime is assumed to be: (a) 1.2x10^'*, (b) 2.5x10'^'*, (c) 5.0x10'^*, 
and (d) 1.0x10^^ years. Dots show the combined signal plus background events and the hatched histograms 
are atmospheric neutrino background events. 

two-body decays. Furthermore, when a proton in an oxygen nucleus decays, the proton hole is 
filled by de-excitation of another proton, resulting in 7 ray emission. The probability of a 6 MeV 7 
ray being emitted is about 40%. This 6 MeV 7 is a characteristic signal used to identify a proton 
decay and to reduce the atmospheric neutrino background. There are three established methods 
for the p —7- T'K^ mode search [12j: (1) look for single muon events with a de-excitation 7 ray just 
before the time of the muon, since the 7 ray is emitted at the time of K'^ production; (2) search 
for an excess of muon events with a momentum of 236 MeV/c in the momentum distribution; and 
(3) search for vr'^ events with a momentum of 205 MeV/c. 

In method (1), the p — )• T'K^{+'y*), — )• + ly candidate events are selected with the 
following criteria: (B-1) a fully contained (FC) event with one ring, (B-2) the ring is ;U-like, (B-3) 
there is a Michel decay electron, (B-4) the muon momentum is between 215 and 260 MeV/c, (B-5) 
the distance of the vertices between the muon and the Michel electron is less than 200 cm, (B-6) 
the time difference between the 7 and the muon is less than 75 ns (~ 6r^+), and (B-7) the number 
of PMTs hit by the 7 is between 4 and 30 to (select 6 MeV energy) . The prompt 7 hits are searched 
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FIG. 43. The significance of the p — > e+Tr" signal derived from counting the events in the signal region 



(Fig. 42) following a 5.6 Megaton-year exposure, as a function of the true proton lifetime. The upper line 



indicates 5 a significance and the lower line indicates 3 a. 

for by sliding a 12 ns wide timing window just before the time of the muon signal as shown in 



Fig. 45 



In method (2), the previous selection criteria are relaxed; the events are not to satisfy (B-5)- 
(B-7), but rather only (B-l)-(B-4). An excess of muon signals in the momentum distribution is 
then searched for by fitting the data with the proton decay signal expectation over the atmospheric 
neutrino background events. In method (3), the vr^ events with a momentum of 205 MeV/c are 
selected. Although the tt"*" does not make a clear Cherenkov ring due to its low momentum, hit 
activity caused by vr^ in the opposite direction of the VfO is used to identify the — t- vr^vr 
signal. The following selection criteria are used: (C-1) FC events with two rings (from vr*^ — t- 77), 
(C-2) both rings are e-like, (C-3) one Michel decay electron from the muon produced by 7r+ — )• /x"*", 
(C-4) the invariant mass of two rings is between 85 and 185 MeV/c^ (vr*^ mass region), (C-5) the 
reconstructed vr'^ momentum is between 175 and 250 MeV/c, (C-6) the visible energy opposite the 
vr*^ (140 - 180 degrees) is between 7 and 17 MeV, and (C-7) the visible energy 90 - 140 degrees from 
vr*^ direction is less than 12 MeV. 

The detection efficiencies are calculated to be 7.1% for method (1), 43% for method (2), and 
6.7% for method (3). The background rates from atmospheric neutrinos are 1.6, 1940, and 6.7 



events/Megaton-year for the methods (1), (2), and (3), respectively. Fig. 46 shows the number 
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FIG. 44. Sensitivity of the proton decay search for the p e+7r" mode as a function of detector exposure. 
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FIG. 45. Schematic view of the prompt 7 hit search in the p — >■ (+7* ) , i4r+ — >• /z"*" + mode. 



of hits distributions in the prompt 7 hits search with all other cuts (B-l)-(B-6) of method (1) 
applied. The running time of Hyper-Kamiokande is assumed to be 10 years which corresponds 
to a 5.6 Megaton- year exposure. The proton lifetime is taken to be: (a) 3.9x10^^ years (current 
limit), (b) 6.3x10^^ years, (c) 1.0x10^^ years, and (d) 1.5x10^^ years. The number of atmospheric 
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neutrino background events is estimated to be 9.0 events. With an assumed background uncertainty 
of 33% (the same as Super-K), proton decay signals from the p — i'K'^{+j*), — ^+ + mode 
could be separated from the background events with a 5 u significance for a proton lifetime of 

U). 



6.3x10^'^ years as shown in Fig. 46 'b), and with 3 a for 1.0x10^^ years as shown in Fig. 



46 



Fig. 47 shows the 90% CL sensitivity curve for the p — t- vK^ mode, by combining all three methods, 
as a function of the detector exposure. In this case, a 1.2 Megaton-year exposure is needed to reach 
a lifetime sensitivity of 10^^ years. That corresponds to two years running Hyper-Kamiokande, or 
53 years of livetime for Super-Kamiokande. 




FIG. 46. Number of hits distributions in the prompt 7 hits search for MC events passing cuts (B-l)-(B-6) 
for a 5.6 Megaton-year exposure. The proton lifetime is assumed to be: (a) 3.9x10'^^, (b) 6.3x10'^'^, (c) 
1.0x10^4, and (d) l.SxlO^-* years. Dots show the signal plus background events and the hatched histograms 
are the atmospheric neutrino background events. 



3. Sensitivity study for other nucleon decay modes 



Although the p — )• e^vr mode is predicted to be the dominant decay mode in many GUT 
models, the branching ratios of other nucleon decay modes are not small. Table XIII shows several 
proton decay modes with the branching ratios predicted by several GUT models. The ratio of 
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FIG. 47. Sensitivity plot for the p vK^ mode as a function of detector exposure. 

neutron to proton lifetimes is also shown in the table. Since we do not know which model (and 
there are many others) is correct, in order to discover proton decay - as well as to examine and 
experimentally constrain the various GUT models - searching for a variety of nucleon decay modes 
is clearly important. 

TABLE XIII. Branching ratios for proton decay and the ratio of neutron to proton lifetimes predicted by 
the SU(5) and SO(IO) models. . 





Br.(%) 






SU(5) 




SO(IO) 


References 


m 


m 


m 


m 


m 


p — > e+TT*^ 


33 


37 


9 


35 


30 


p — > e+ry" 


12 


7 


3 


15 


13 


p — > e+/3° 


17 


2 


21 


2 


2 


p — > e+cj*^ 


22 


18 


56 


17 


14 


Others 


17 


35 


11 


31 


31 




0.8 


1.0 


1.3 







The Hyper-Kamiokande sensitivities for other nucleon decay modes are estimated based on the 
efficiencies and background rates of SK-II . Table XIV shows the 90 % CL sensitivities with a 
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5.6 Megaton-year exposure (10 years running of Hyper-Kamiokande) . The current lifetime limits 
are also shown in the table for reference. In all cases, we could explore an order of magnitude 
longer lifetime regions. 

The decay modes in Table |XIV| all conserve the baryon number minus the lepton number, 
{B — L). Recently another {B + L) conserving mode, n — t- e~K^ , was also given attention and 
searched for by Super-Kamiokande. In n — )• e~K~^, the stops in the water and decays into 

+ v. The final state particles observed in n — )• e~K~^, — )• (j.'^v are e~ and /i"*". Both e~ and 

have monochromatic momenta as a result of originating from two-body decays. Furthermore, 
the timing of /i^ events are delayed with respect to the e~ events because of the lifetime. In 
SK-II, the estimated efficiencies and the background rate are 8.4% and 1.1 events/Megaton-year, 
respectively. From those numbers, the sensitivity to the n — ?■ mode with a 5.6 Megaton-year 

exposure is estimated to be 2.2x10^"^ years. 

The possibility of nn oscillation is another interesting phenomenon; it violates the baryon 
number [B) by |Ai?| = 2. These nn oscillations have been searched for in Super-Kamiokande 
with a 0.09 Megaton-year worth of data [76j. Further improvement of the nn oscillation search is 
expected in Hyper-Kamiokande. 

TABLE XIV. Summary of the sensitivities from a 5.6 Megaton-year exposure compared with current life- 
time limits. The current limits for p — > g+tt*^, p — > /^^tt^, and p — > VK^ arc the result of analyzing a 
0.22 Megaton-year data sample, while the other modes are from a 0.14 Megaton-year exposure. 



Mode 


Sensitivity (90% CL) 


Current limit 


p — > e+TT*^ 


13x10^4 years 


1.3x1034 years 


p /x+7r° 


9.0x1034 


1.1x1034 


p — > e+ry" 


5.0x10^4 


0.42x1034 


p — > ii^rf^ 


3.0x10^'' 


0.13x1034 


p — > 


1.0x1034 


0.07x1034 


p — > 


0.37x1034 


0.02x1034 


p — > e'^LO^ 


0.84x1034 


0.03x1034 


p — > ^.'^uP 


0.88x1034 


0.08x1034 


n — > e+TT^ 


3.8x1034 


0.20x1034 


n — > /i+TT^ 


2.9x1034 


0.10x1034 


p — >• VK+ 


2.5x1034 


0.40x1034 
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4-. Summary 

We have studied the sensitivity of various nucleon decay searches in Hyper-Kamiokande. Ta- 



ble XV shows the summary of the study for the highhghted modes, p — )• e'^ir^ and p — )• uK~^. If the 
proton Ufetime is shorter than 5.7x10^'^ years for the p — >• e'^n^ mode, or shorter than 1.0x10^^ 
years for p — )• uK~^ , we could discover a signal over the atmospheric neutrino background events 
with a 3a significance by collecting a 5.6 Megaton-year exposure. 

TABLE XV. Summary of the sensitivity study for a 5.6 Megaton-year exposure for the p — !■ e^Tr" and 
p VK+ modes. For p VK+ , method (1) n + 6MeV7 is labeled "Meth.l", method (2) (/x) as "Meth.2", 
and method (3) tt+tt^ as "Meth.3". 





p — > e+Tr" 


p vK+ 






Meth.l 


Meth.2 


Meth.3 


Efficiency (%) 


45 


7.1 


43 


6.7 


Background (/Mton-yr) 


1.6 


1.6 


1940 


6.7 


90% Sensitivity (xlO^* yrs) 


13 


2.5 


3(7 Discovery potential ( x 10'^'* yrs) 


5.7 


1.0 



Fig. 48 compares the proton lifetime predictions made by several leading theoretical models with 
the current experimental limits from Super-Kamiokande and the sensitivity of a 10 year Hyper- 
Kamiokande run. Although Super-Kamiokande has been running well and gives us the most 
stringent limits for various nucleon decay searches, we have not observed any evidence of nucleon 
decay yet. The sensitivity of Hyper-Kamiokande can cover most of the predicted range of the 



major GUTs models. Fig. 49 shows the sensitivity of Super-Kamiokande and Hyper-Kamiokande 
as a function of the year. Hyper-Kamiokande is assumed to start data-taking from 2019. The 
nucleon decay results of Hyper-Kamiokande can overtake Super-Kamiokande within a single year 
of running. There is no question that to explore an order of magnitude longer lifetime regions, a 
larger detector is absolutely necessary. As seen in this sensitivity study, Hyper-Kamiokande will 
open up a new decade, in more ways than one, in the search for nucleon decay. 
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FIG. 48. Proton lifetime predictions of several GUT models, the current experimental limits (90% CL) by 
Super-K, and the sensitivities of Hyper-Kamiokande with a 5.6 Megaton-year exposure. Hyper-Kamiokande 
can cover most of the predicted range of the leading GUT models. 




FIG. 49. The proton decay search sensitivity as a function of year. The left plot is for the p — > e+Tr" mode 
and the right is for the p — VK^ mode. Hyper-Kamiokande is assumed to start from 2019; its results will 
overtake the Super-Kamiokande limits within one year. 
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D. Solar neutrinos 



The effective fiducial volume for the solar neutrino analysis in Hyper-Kamiokande will be in- 
creased by a factor of 27 compared to Super-Kamiokande. In this section, high statistics measure- 
ments of solar neutrinos in Hyper-K are discussed. 

In solar neutrino oscillations, regeneration of the electron neutrinos through the Mikheyev- 
Smirnov-Wolfenstein (MSW) matter effect |65H67) in the Earth is expected. Regeneration of the 
solar electron neutrinos in the Earth would constitute concrete evidence of the MSW matter effect, 
and so it is important to experimentally observe this phenomenon. However, the matter effect has 
not been directly confirmed yet, since the sensitivities of the current solar neutrino experiments 
are not sufficient. According to the MSW model, the observed solar neutrino event rate in water 
Cherenkov detectors in the nighttime is expected to be higher - by about a few percent in the 
current solar neutrino oscillation parameter region - than that in the daytime. We would like to 
measure this difference in Hyper-Kamiokande. 

Hyper-K also could be used for variability analyses of the Sun. For example, the solar 
neutrino flux highly depends on the Sun's present core temperature. Unlike multiple scattered, 
random-walking photons or slow-moving helioseismic waves, free streaming solar neutrinos are the 
only available messengers with which to precisely investigate ongoing conditions in the core region 
of the Sun. Hyper-K, with its unprecedented statistical power, could measure the solar neutrino 
flux over short time periods. Therefore, short time variability of the temperature in the solar core 
could be monitored by the solar neutrinos in Hyper-K. 

In order to achieve these precision measurements, background event levels must be sufficiently 
small. Here, we have estimated the basic performance of Hyper-K for low energy events assuming 
some typical background levels. In this study, the current analysis tools and the detector simulation 
for the low energy analysis [16] in Super-K were used. The dark rate of the PMTs and the water 
transparency were assumed to be similar to those in the current Super-K detector. A brief summary 



of the low energy event reconstruction performance in Hyper-K is listed in Table IV 



The analysis threshold of the total energy of the recoil electrons in Hyper-K will be 7.0 MeV or 
lower, since a 7.0 MeV threshold was previously achieved in the SK-II solar neutrino analysis |77] . 
The current analysis tools will work all the way down to 4.5 MeV in Hyper-K with a vertex 
resolution of 3.0 m. Not surprisingly, higher energy events will be reconstructed with even better 
vertex resolution. 
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1. Background estimation 

The major background sources for the solar neutrino measurements are the radioactive 
daughter isotopes of Rn-222 in water and the radioactive spallation products created by cosmic- 
ray muons. Rn-222 will be reduced to a similar (or lower) level as that currently in the Super-K 
detector, since Hyper-K will employ a similar water purification system and design improvements 
may well occur over the next several years. However, the spallation products will be increased in 
Hyper-K. In the current design, the cosmic-ray muon rate is expected to be increased by a factor 



of about 10 in equal volumes, as discussed in Sec. HA The spallation products will not simply be 



increased by the same factor. This is because high energy cosmic-ray muons tend to produce the 
spallation products, while the average energy of the cosmic-ray muons at the shallower Hyper-K 
site is expected to be lower than that at the deeper Super-K site; greater overburden means less 
muons, but it also means those that do get through are more energetic. We have estimated the 
average energies of the cosmic-ray muons to be ~ 560 GeV at the Super-K site and ~ 300 GeV at 
the Hyper-K site. Considering that the spallation production cross section is proportional to the 
0.7-th power of the cosmic-ray muons' energy, the density of spallation products will be increased 
by a factor of 6-7. We found the remaining spallation products will be increased by another factor 
of 3 (at most) with the current analysis tools. This is due to decreasing efficiency for separating 
spallation products from signal events with increasing the cosmic-ray muon rate. So, the density of 
the remaining spallation products will be increased by a factor of 20 at most in Hyper-K. However, 
this could (and most likely will) be reduced by ongoing improvements of the analysis tools. 

In Super-K, angular information is used to extract the solar neutrino signal events [78j. We 
have estimated the possible effect of the background level in the signal extraction after considering 
angular information. 

In this study, we used 9.0-9.5 MeV Super-K-I data as a reference. The extracted solar 
neutrino signal events and background events in this energy region over the entire run period 
(0.09 Megaton-years) were 1350 events and 7700 events, respectively. So, the Signal-to-Noise 
(S/N) ratio is 18%. We made artificial data samples with reduced S/N ratios, then applied the 
signal extraction. As a result, we found the expected statistical error is almost the square root 
of 2 ~ 15 times the number of signal events for 1 ~ 20 times the Super-K-I background level, 



respectively. Table XVI shows a summary of the expected statistical errors in a Super-K-I type 
detector with increased backgrounds, as well as that of Hyper-K. Once the angular distribution is 
used to extract the solar signal, the statistical error on this signal would be reduced by a factor 
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TABLE XVI. Expected statistical uncertainties for 10000 signal events with increased background levels. 
The Super-K-I solar neutrino data sample between 9.0-9.5 MeV was used as a reference The 3rd column is 
the Hyper-K factor relative to Super-K given the same observation time. To estimate the 3rd column, the 
same detector resolution and 0.56 Mton fiducial volume are assumed in Hyper-K. 

Background level Stat. err. in SK Stat. err. in HK 



SKTBGx20 3.6% x 1/2.0 

SKTBGxlO 2.7% x 1/2.7 

SKT BG x7 2.4% x 1/3.1 

SKT BG x5 2.1% xl/3.5 

SKTBG 1.4% X 1/5.2 



of 2.0 in Hyper-K, even though the background level is increased by a factor of 20, for the same 
observation time assuming both detectors have identical resolution. 

In summary, Hyper-K will provide higher statistical measurements of solar neutrinos than Super- 
K, even though there will be more spallation backgrounds. 



2. Oscillation study 

In solar neutrino oscillations, a difference in the solar neutrino event rates during the daytime 
and the nighttime is expected from the MSW effect in the Earth. This is called the day/night 
asymmetry; it has not yet been observed. In Hyper-K, a precise measurement of the day/night 
asymmetry will be performed using higher statistics than those available in Super-K. 

The upper plots in Fig. |50| show the expected day/night asymmetries with different lower energy 
thresholds. The expected day/night asymmetry is at about the 1% level around the current solar 
global oscillation parameters. In order to observe the day /night asymmetry in Hyper-K, we must 
reduce the up-down systematic uncertainty below that level. 

The expected day/night asymmetry in the high energy region is larger than that in the low 



energy region, as shown in the lower left plot in Fig. 50 So, high statistics data in this higher energy 
region would be desirable. We have studied two typical values of systematic uncertainty, where the 
one at 1.3% corresponds to the up-down systematic uncertainty of Super-K. From the lower right 



plot in Fig. 50, the most sensitive lower energy threshold would be 6 MeV and 8 MeV for the 0.5% 



and 1.3% up-down systematic uncertainties, respectively. Figure 51 shows expected day/night 
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(N-D)/(N+D)*2 5.0-20MeV (N-D)/(N+D)*2 8.0-20MeV 




Energy threshold (MeV) Energy threshold (MeV) 

FIG. 50. Expected day/night asymmetry in a megaton water Cherenkov detector. 40% photo-coverage, 
0.5 Megaton-years daytime data and 0.5 Megaton-years nighttime data are assumed. The effect of back- 
ground events and reduction efficiencies are not considered. Upper left: expected day/night asymmetry 
in the 5.0-20 MeV electron total energy region. Upper right: expected day/night asymmetry in the 8.0- 
20 MeV region. Lower left: expected day/night asymmetry with uncertainties as a function of the lower 
energy threshold at (tan^ 6*12, Am^ J = (0.40,7.9 x lO^^eV^). The upper energy threshold is 20 MeV. 
The meaning of the different colors are defined in the lower right plot. Lower right: expected day /night 
significance as a function of the energy threshold. 

significance as a function of the observation time. Since the expected day/night asymmetry is 
small, it will be important to reduce the systematic uncertainties in order to observe the day/night 
asymmetry with high precision. We believe that this should be possible, especially if we design 
and prepare the necessary calibration devices during detector construction. 

3. Time variation study 

Solar neutrinos could be used as a direct probe of the nuclear reactions taking place in the 
solar core. In particular, the *B solar neutrino flux has a remarkable T^^ dependence according 
to Standard Solar Model (SSM) [79]. Here, T is the solar core temperature, and with such a 
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FIG. 51. Expected day/night significance as a function of the observation time near the solar global os- 
cillation best-fit parameters. The Super-K-I value of S/N is assumed. The total electron energy region is 
5.0-20 MeV. 

high-order dependence it is possible that even modest changes in the solar core temperature could 
be ampHfied into something detectable via measurements of the solar neutrino flux. 



Assuming the statistical uncertainties estimated in Sec. Ill D 1 can be used for Hyper-K, the 
expected uncertainty on the solar core temperature when the background level is increased by a 
factor of 20 would be the following: 

ctt 1 Cat 

~T ~ \^~W ~ IS-N 

Here N, a^, and ctn are the number of observed '^B solar neutrinos, error in T, and error in A'", 
respectively. The expected number of observed *B solar neutrinos in Hyper-K is 200 events per 
day above 7.0 MeV, as shown in Table [lT| When is 200, cjt/T will be 0.015. Therefore, the solar 
core temperature could be monitored within a few percent accuracy day by day. Naturally, by 
integrating over longer periods, more subtle temperature changes - potentially down to the 0.1% 
level - could be monitored. 



4-. Summary 

In this section, rough estimates of potential solar neutrino measurements are reported. The 
solar neutrino analysis is sensitive to the detector resolutions and background levels. We have 
estimated expected sensitivities based on the current Super-K analysis tools. 
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As a result of its shallower site, the increase of the background level in Hyper-K will be up to 
a factor of 20 as compared to Super-K. However - due to its much greater size - the statistical 
uncertainties on solar neutrino measurements would actually be reduced by a factor of at least two 
in Hypcr-K as compared to Super-K on an equal time basis, assuming similar detector resolutions. 

The day /night asymmetry of the solar neutrino flux - concrete evidence of the matter effect on 
oscillations - could be discovered and then precisely measured in Hypcr-K, given that the detector 
up-down response is understood to better than about 1%. Good calibration tools will be a must 
for this physics. 

Hyper-K will provide short time and high precision variability analyses of the solar core activity. 
The solar core temperature could be monitored within a few percent accuracy day by day, and to 
a tenth of a percent over the period of several months. 
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E. Astrophysics 

In this section, potential signals arising from several astrophysical sources are considered, and 
the prospects for their detection and study in Hyper-Kamiokande are briefly discussed. 

1. Supernova burst neutrinos 

Core collapse supernova explosions are the last process in the evolution of massive (> SMgun) 
stars. Working their way successively through periods of predominantly hydrogen fusion, helium 
fusion, and so on, eventually silicon fusion starts making iron. Once an iron core has formed, 
no more energy can be released via its fusion into still-heavier elements, and the hydrodynamic 
balance between gravity and stellar burning is finally and catastrophically disrupted. The sudden 
gravitational collapse of their iron cores - each one of which then goes on to form either a neutron 
star or a black hole - is the main source of energy from this type of supernova explosion. The energy 
released by a supernova is estimated to be '-^ 3 x 10^^ ergs, making it one of the most energetic 
phenomena in the universe. Since neutrinos interact weakly with matter, almost 99% of the released 
energy from the exploding star is carried out by neutrinos. As a result, the detection of supernova 
neutrinos gives direct information of energy flow during the explosion. The neutrino emission from 
a core collapse supernova starts with a short (~10 millisecond) burst phase of electron captures 
{p + e~ ^ n + Ue) called the neutronization burst, which releases about 10^-^ ergs. Following that, 
the majority of the burst energy is released by an accretion phase (<~'l second) and a cooling 
phase (several seconds) in which all three types of neutrinos (including anti-neutrinos) are emitted. 

The observation of a handful (25 in total) of supernova burst neutrinos from SN1987a by the 
Kamiokande, 1MB, and Baksan experiments proved the basic scenario of the supernova explosion 
was correct. However, more than two decades later the detailed mechanism of explosions is still not 
known. Several groups around the world are trying to explode supernovae in computer simulations. 
However, they have not been completely successful yet. It seems that some physical processes 
arc lacking in the present simulations. In order to fully understand supernova explosions, it is 
necessary to detect many more supernova neutrinos. If a supernova explosion occurs halfway 
across our galaxy, the Hyper-Kamiokande detector would detect approximately 170,000^^260,000 
neutrino events. This very large statistical sample should at last reveal the detailed mechanism of 
supernova explosions. 

Supernova neutrinos also give us an opportunity to investigate unknown properties of neutrinos. 
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Because of quite high densities at the supernova core, neutrino-neutrino self-interaction is not 
neghgible. As described later in this section, the neutrino mass hierarchy could be determined 
using the time variation of the energy spectrum. The very short rise time of the outbreak of the 
burst enables us to discuss neutrino direct masses. Also, a detection of neutrinos from distant 
sources give us a chance to check the lifetime of neutrinos. Indeed, by far the oldest neutrinos ever 
seen are those of SN1987a. They were about 170,000 years old when they arrived, while the next 
oldest would be the solar neutrinos, seen on Earth a mere eight minutes after their birth in the 
heart of the Sun. 

Taking into account the Mikheyev-Smirnov-Wolfenstein (MSW) matter effect through the stellar 
medium, the flux of each neutrino type emitted from a supernova is related to the originally 
produced fluxes {Fj^^, F^^ and F^^, where Vx is t'/^.r and i'^^t) by the following formulas [801 181] • 
For normal hierarchy, 

F^, ~cos2^i2F^°, +sin2ei2< , 

F,^ ~ sin^ ^i2PhF° + (1 - sin^ ^laP/^X, 
F., + F,^ ~ (1 - sin2 ei2PH)Fl + (1 + sin2 euPH)F!^^, 
Fv, + Ft,^ ^{l- cos^ ei2)Fl + (1 + cos2 ei2)Fl, 

and, for inverted hierarchy, 

~ cos^ ei2PHFl + (1 - cos^ ei2PH)Fl , 
F,, ~sin2 0i2F° +cos2 0i2<, 
F,^ + F,^ :^ (1 - sin2 ^12)^° + (1 + sin^ ^12)^°^, 
Fu, + F^. =^ (1 - cos2 ei2PH)Fl + (1 + sin2 0i2Fh)FO^, 

where Ph is the crossing probability through the matter resonant layer corresponding to 
Ph = ^{Ph = 1) for adiabatic (non-adiabatic) transition. For sin^ 2^13 > 10^^, adiabatic transition 
is expected for the matter transition in the supernova envelope. The supernova neutrino spectrum 
is affected not only by stellar matter but also by other neutrinos and anti-neutrinos at the high 
density core (so-called collective effects). The collective effects swap the Ve and z/e spectra with 
those of i^x in certain energy intervals bounded by sharp spectral splits [82] . The combination of the 
collective effects and the stellar matter effects might weaken swapping the energy spectra of fe/z^e 
and Ux- So, in the following description of the performance of the Hyper-K detector, three cases 
are considered in order to fully cover the possible variation of expectations: (1) no oscillations, (2) 
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FIG. 52. Expected time profile of a supernova at 10 kpc. Left, center, and right figures show profiles for no 
oscillation, normal hierarchy (N.H.), and inverted hierarchy (I.H.), respectively. Black, red, purple, and light 
blue curves show event rates for interactions of inverse beta [ve^V ^ e^"+n), j/e-scattering(i^+e^ — > j/+e^), 
cC(i/e + ^ e- + i6f(*)), and ^^Q) CC [y^ + ^'^O ^ e+ + respectively. The numbers 

in parentheses are integrated number of events over the burst. The fluxes and energy spectrums are from 
the Livermore simulation 1831 



normal hierarchy with Pn = 0, and (3) inverted hierarchy with Pu = 0. An assumed oscillation 
parameter relevant for the description is sin^ 9i2 = 0.31. Concerning the neutrino fluxes and energy 
spectrums at the production site, we used results obtained by the Livermore simulation |83j . 



Figure 52 shows time profiles for various interactions expected at the Hyper-K detector, if a 
supernova at a distance of 10 kiloparsecs (kpc) is observed. This distance is a bit farther than the 
center of the Milky Way galaxy at 8.5 kpc; it is chosen as being representative of what we might 
expect since a volume with a radius of 10 kpc centered at Earth includes about half the stars in 
the galaxy. The three graphs in the figure show the cases of no oscillation, normal hierarchy (N.H.) 
and inverted hierarchy (I.H.), respectively. Colored curves in the figure show event rates for inverse 
beta (Pe +P e+ + n), z^e-scattering(i/ + e~ ^ u + e~), Ve ^^O CC{ve + ^'^O e~ + ^^pC*)), 
and Ve ^''O CC {ue + ^^O e+ + 16n(*)). The burst time period is about 10 s and the peak 
event rate of inverse beta events reaches about 200 kHz. A sharp timing spike is expected for 
z^e-scattering events at the time of neutronization, which will be discussed in more detail later. 



Visible energy distributions of each interaction are shown in Fig. 53 where the visible energy is 
the electron-equivalent energy measured by a Cherenkov detector. The distribution of inverse beta 
events directly gives the energy spectrum of i>(. because the energy of the positrons (-E'e+ ) is simply 
— 1.3 MeV. The energy spectrum of i/e-scattering events has enhancement at lower energies 
because of the nature of the interaction. 
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visible energy (MeV) 

FIG. 53. Visible energy spectrum for each interaction for a supernova at 10 kpc. Black, red, purple, and 
light blue curves show event rates for interactions of inverse beta, i^e-scattering, I'l^O CC, and D^^O CC, 
respectively. Solid, dashed, and dotted curves correspond to no oscillation, N.H., and I.H., respectively. 



Figure 54 shows the expected number of supernova neutrino events at Hyper-K versus the 
distance to a supernova. At the Hyper-K detector, we expect to see about 165,000~230,000 inverse 
beta events, 7,000~8,000 z^e-scattering events, 300~14,000 ufO CC events, and 2,000~13,000 i?l^O 
CC events, in total 170,000~260,000 events, for a 10 kpc supernova. The range of each of these 
numbers covers possible variations due to the neutrino oscillation scenario (no oscillation, N.H., 
or I.H.). Even for a supernova at M31 (Andromeda Galaxy), about 30~50 events are expected at 
Hyper-K. In the case of the Large Magellanic Cloud (LMC) where SN1987a was located, about 
7,000~10,000 events are expected. 



Figure 55 shows expected angular distributions with respect to the direction of the supernova 
for four visible energy ranges. The inverse beta events have a nearly isotropic angular distribution. 
On the other hand, z^e-scattering events have a strong peak in the direction coming from the 
supernova. Since the visible energy of z^e-scattering events are lower than the inverse beta events, 
the angular distributions for lower energy events show more enhanced peaks. The direction of a 
supernova at 10 kpc can be reconstructed with an accuracy of about two degrees according to these 
angular distributions. In Hyper-K, we can statistically extract an energy distribution of z^e + t^x 
{X = /i, r) events using the angular distributions in much the same way as solar neutrino signals 
are separated from background in Super-K. Although the effect of neutrino oscillations must be 
taken into account, the + vx spectrum gives another handle for discussing the temperature 
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FIG. 54. Expected number of supernova burst events for each interaction as a function of the distance to a 
supernova. The band of each hue shows the possible variation due to the assumption of neutrino oscillations. 
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FIG. 56. Inverse beta event rate (left) and mean energy of D,, (right) predicted by supernova simulations 

[831 - 187] for the first 0.3 seconds after the onset of a 10 kpc distant burst. 



of neutrinos. Note that inverse beta events directly provide a very precise measurement of the 
temperature of I'e- Hyper-K will be able to evaluate the temperature difference between and 
i^e + t^x- This would be a valuable input to model builders. 



Figure 56 shows inverse beta event rates and mean Uf. energy distributions predicted by various 
models |83H87j for the first 0.3 sec after the onset of a burst. The statistical error of Hyper-K is 
much smaller than the difference between the models, and so Hyper-K should give crucial data 



for comparing model predictions. The left plot in Fig. 56 shows that about 300-1000 events are 
expected in the first 20 millisecond bin. This means that the onset time can be determined with 
an accuracy of about 0.03 ms. This is precise enough to allow examination of the infall of the core 
in conjunction with the signals of neutronization (see below) as well as possible data from future 
gravitational wave detectors. 

We can also use the sharp rise of the burst to make a measurement of the absolute mass of 
neutrinos. Because of the finite mass of neutrinos, their arrival times will depend on their energies. 
This relation is expressed as 



At = 5.15 msec 



m \ ' 
Yev) 



10 MeV ; 



(7) 



,10 kpc, 

where At is the time delay with respect to that assuming zero neutrino mass, D is the distance to 
the supernova, m is the absolute mass of a neutrino, and Ei, is the neutrino energy. Totani |88| 
discussed Super-Kamiokande's sensitivity to neutrino mass using the energy dependence of the 
rise time; scaling these results to the much larger statistics provided by Hyper-K, we expect a 
sensitivity of 0.5~1.3 eV for the absolute neutrino mass [89] • Note that this measurement of the 
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FIG. 57. Expected event rate at the time of neutronization burst for a supernova at 10 kpc. Red and green 
show event rates for i/e-scattering and inverse beta events, respectively. Sohd, dotted, and dashed curved 
indicate the neutrino osciUation scenarios of no osciUation, N.H., and I.H., respectively. 



absolute neutrino mass does not depend on whether the neutrino is a Dirac or Majorana particle. 
A sharp timing peak at the moment of neutronization is expected in z/e-scattering events as 



shown in Fig. 57, The expected number of z^e-scattering events at the neutronization burst is ~20, 
~56, and ~130 for N.H., I.H., and no oscillation, respectively, for a supernova at 10 kpc. Although 
the number of inverse beta events is ~345 (N.H.), ~700 (I.H.), and ~190(no oscillation) in the 
10 ms bin of the neutronization burst, the number of events in the direction of the supernova is 
typically 1/10 of the total events. So, the ratio of signal events (z^e-scattering) to other events 
(inverse beta) is expected to be about 20/33 (N.H.), 52/70 (I.H.) and 130/19 (no oscillation). 
Thus, the ve scattering events can be identified with high statistical significance thanks to the 
directionality of z^e-scattering. 

Neutrino oscillations could be studied using supernova neutrino events. There are many papers 
which discuss the possibility of extracting signatures of neutrino oscillations free from uncertainties 
of supernova models (HU [90ti98] . One big advantage of supernova neutrinos over other neutrino 
sources (solar, atmospheric, accelerator neutrinos) is that they inevitably pass through very high 
density matter on their way to the detector. This gives a sizable effect in the time variation of 
the energy spectrum even for small sin^ ^13 [81^ [91] I99j . As an example, figures from the paper 



by Fogli et al. [SI] are shown in Fig. 58, The propagation of the supernova shock wave causes 
time variations in the matter density profile through which the neutrinos must travel. Because of 
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FIG. 58. Time variation of the neutrino event rate affected by neutrino conversion by matter due to shock 
wave propagation (reproduced from 81 ). Left (right) plot shows the time variation of inverse beta events 
for the energy range of 20±5 MeV (45±5 MeV). Sohd black, dashed red, and blue dotted histograms show 
the event rates for I.H. with shock wave propagation, I.H. with static matter density profile, and N.H., 
respectively. It has been assumed that sin^ 6*13 = 10^^. 



neutrino conversion by matter, there might be a bump in the time variation of the inverse beta 



event rate for a particular energy range (i.e., 45ib5 MeV as shown in Fig. 58 right)) while no change 
is observed in the event rate near the spectrum peak (i.e., 20ib5 MeV as shown in Fig. [58|left)). 
This effect is observed only in the case of inverted mass hierarchy; this is one way in which the 
mass hierarchy could be determined by a supernova burst. 

In Hyper-K, it could be possible to detect burst neutrinos from supernovae in nearby galaxies. 
As described above, we expect to observe a very large number of neutrino events from a galactic 
supernova. However, galactic supernovae are expected to happen once per 30-50 years. So, we 
cannot count on seeing many galactic supernova bursts. In order to examine a variety of supernova 
bursts, supernovae from nearby galaxies are useful even though the expected number of detected 
events from any single such burst are small. Furthermore, in order to fully understand the spectrum 
of supernova relic neutrinos (see next sub-section), collecting an energy spectrum without the 
complications of varying red-shift effects is highly desirable. The supernova events from nearby 
galaxies provide a reference energy spectrum for this purpose. The supernovae in nearby galaxies 



was discussed by S. Ando, J. F. Beacom and H. Yiiksel [TOO]; a figure from their paper is shown 



in Fig. 59 left). It shows the cumulative supernova rate versus distance and indicates that if 
Hyper-Kamiokande can see signals out to 4 Mpc then we could expect a supernova about every 
three years. It should also be noted that the paper says recent, more sensitive astronomical 
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Distance D [Mpc] Distance (Mpc) 



FIG. 59. (Left) Cumulative calculated supernova rate versus distance for supernovae in nearby galaxies 
(reproduced from jlOQj). (Right) Detection probability of supernova neutrinos versus distance at Hyper-K 
assuming a 0.56 Megaton fiducial volume and 10 MeV threshold for this analysis. Black, green, and blue 
curves show the detection efficiency resulting in requiring at least or equal to one, two, and three events per 
burst, respectively. Solid, dotted, and dashed curves are for neutrino oscillation scenarios of no oscillation, 
N.H., and I.H., respectively. 



observations indicate that the true nearby supernova rates are probably about 3 times higher than 



this conservative calculation. Figure 59 fright) shows detection probability versus distance for the 
Hyper-K detector. In this estimate, energy of neutrino events is required to be more than 10 MeV 
and the vertex position of the events should be within the fiducial volume (0.56 Megatons). If we 
require the number of events to be more than or equal to one(two), the detection probability is 
52~69%(17~33%) for a supernova at 4 Mpc. If we can use a tight timing coincidence with other 
types of supernova sensors (e.g. gravitational wave detectors), we should be able to identify even 
single supernova neutrinos. Assuming the observed supernova rate in nearby galaxies, we expect 
to collect about 10-20 supernova neutrino events from them during 20 years of operating Hyper-K. 



2. Supernova relic neutrinos 

There are about 10^" stars in the universe (~10^'' galaxies in the universe, and each galaxy 
has about 10^*^ stars). Because about 0.3% of the stars have masses larger than 8 times the solar 
mass, it is estimated that 10^^ supernova explosions have occurred over the entire history of the 
universe. This means that on average one supernova explosion has been occurring every second 
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FIG. 60. Predictions of the supernova relic neutrino (SRN) spectrum. Fluxes of reactor neutrinos and 
atmospheric neutrinos are also shown. 



somewhere in the universe. The neutrinos produced by all of the supernova explosions since the 
beginning of the universe are called supernova relic neutrinos (SRN). They must fill the present 
universe and their flux is estimated to be a few tens/cm^/sec. If we can detect these neutrinos, it 
is possible to explore the history of how heavy elements have been synthesized since the onset of 



stellar formation. Figure 60 shows the SRN spectra predicted by various models. Although searches 
for SRN have been conducted at large underground detectors, no evidence of SRN signals has yet 
been obtained because of the small flux of SRN. The expected inverse beta {ueP e^n) event 
rate at Super-Kamiokande is 0.8-5 events/year above 10 MeV, but because of the large number 
of spallation products and the low energy atmospheric neutrino background (decay electrons from 
muons below Cherenkov threshold produced by atmospheric muon neutrinos, the so-called invisible 
muon background), SRN signals have not yet been observed at Super-Kamiokande. In order to 
reduce background, lower the energy threshold, individually identify true inverse beta events by 
tagging their neutrons, and thereby positively detect SRN signals at Super-Kamiokande, a project 
to add 0.1% gadolinium to the tank (the GADZOOKS! project) was proposed by J.F.Beacom and 
M.R.Vagins [lOlj : very active R&D work for the project is ongoing. The first observation of the 
SRN could be made by the GADZOOKS! project, but in order to measure the spectrum of the 
SRN and analyze the history of the universe we need a megaton-scale detector. 

Figure |6l] shows expected SRN signals at Hyper-K with 10 years' livetime without tagging 
neutrons. Because of the high background rate below 20 MeV from spallation products, the 
detection of SRN signals is limited to above 20 MeV, while above 30 MeV the atmospheric neutrino 
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FIG. 61. Expected spectrum of SRN signals at Hyper-K with 10 years of livetime without tagging neutrons. 
The black dots show signal+background (red component). Green and blue show background contributions 
from the invisible muon and Ve components of atmospheric neutrinos. Without tagging neutrons, spallation 
background dominates below 20 MeV. The SRN flux prediction by Ando et al. |102| is used, and a signal 
selection efficiency of 90% is assumed. 



backgrounds completely overwhelm the signal. The expected number of SRN events in E = 
20 — 30 MeV is about 310/lOyrs assuming the flux prediction of Ando et al. |102j and an event 
selection efficiency of 90%. The number of background events from atmospheric neutrinos (invisible 
muon and i>e) is estimated to be 2200/lOyrs. So, it is possible to detect SRN signals with high 
statistical significance. 

However, it is important to be able to measure the SRN spectrum down to ~10 MeV in order 
to explore the history of supernova bursts back to the epoch of red shift (z) ~1. Therefore, in the 
following discussion the expected SRN signal with gadolinium neutron tagging is considered. 

Inverse beta reactions can be identified by coincident detection of both positron and delayed 
neutron signals, and requiring tight spatial and temporal correlations between them. With 0.1% by 
mass of gadolinium (Gd) dissolved in the water, neutrons are captured on gadolinium with about 
90% efficiency; the excited Gd nuclei then de-excite by emitting 8 MeV gamma cascades. The 
time correlation of about 20 //sec between the positron and the Gd(n,gammas)Gd cascade signals, 
and the vertex correlation within about 50 cm are strong indicators of a real inverse beta event. 
Requiring both correlations (as well as requiring the prompt event to be Cherenkov-like and the 
delayed event to be isotropic) can be used to reduce background of spallation products by many 
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Vertices within 50cm Energy[MeV] 

FIG. 62. The left plot shows how to tag neutrons produced by the inverse beta reaction, while the right 
plot shows the electron-equivalent energy spectrum of the Gd(n, gammas) Gd gamma cascade as measured 
in the Super-Kamiokande detector |103) . 



orders of magnitude while also reducing invisible muon backgrounds by about a factor of 5. Since 
the 8 MeV gamma cascade of Gd(n, gammas) Gd reaction produces multiple gammas, its visible 



energy is much lower than 8 MeV. This was measured to be 4.3 MeV as shown in Fig. 62 right). 
In order to reconstruct vertex and energy of such events, about 20% photocoverage is necessary. 



Figure 63 shows the expected observation of the SRN spectrum in Hyper-Kamiokande. In this plot, 
an SRN flux prediction by Ando et al. [102] was used, and a detection efficiency of 8 MeV gamma 
cascades of 67% (90% for capture efficiency and 74% for event selection in which vertex coincidence 
of less than 2 m, time correlation less than 60 fisec, energy of delayed signal more than 3 MeV, 
and spherical pattern of the Gd gamma event were required) and a reduction factor of the invisible 
muon background of 5 were assumed. The expected number of SRN events in the energy range 
of 10-30 MeV is about 830 with 10 years of live time. These large statistics are sufficient to start 
discussing the supernova history of the universe. In addition, by comparing of the total SRN flux 
with optical supernova rate observations, a determination of the fraction of failed (optically dark) 
supernova explosions, currently unknown but thought to occur in not less than 5% and perhaps as 
many as 50% of all explosions, will be possible. 

Possible backgrounds to the SRN search are (1) chance coincidences, (2) spallation products 
with accompanying neutrons, and (3) the resolution tail of the reactor neutrinos. For (1) chance 
coincidences, the possible source of the prompt event is the spallation products. By requiring time 
coincidence, vertex correlation and energy and pattern of the delayed event, the chance coincidence 
rate can be reduced below the level of the expected SRN signal by reducing the radioactivity 
of PMTs and their protective cases by a factor of about 5. For (2) spallation products with 
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FIG. 63. Expected spectrum of the SRN signals at Hyper-K with 10 years of hvetime. The black dots show 
signal+background (red component). Green and blue show background contributions from the invisible 
muon and Ve components of atmospheric neutrinos. The SRN flux prediction of Ando et al. ^102j was used, 
and a 67% detection efficiency of 8 MeV gamma cascades and a factor of 5 reduction in the invisible muon 
background were assumed. 



accompanying neutrons, the only possible spallation product is ^Li and an estimation by a GEANT4 

Because of the short half-life of ^Li (ri/2=0.18sec), a high rejection 



simulation is shown in Fig. 



64 



efficiency of ~99.5% is expected. With this expectation, the ^Li background is less than the signal 
level above 12 MeV; this could be lowered by further development of the background reduction 
technique. For (3) the resolution tail of the reactor neutrinos, the estimated background rate is 



about 380(80)/10 years above 10 MeV (11 MeV) as shown in Fig. 64 with full reactor intensity. 
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FIG. 64. Green (red) curve shows the estimated ^Li production rate before (after) applying cuts based 
on a correlation with cosmic ray muons. Blue shows estimated background from reactor neutrinos at full 
intensity. Black data points show expected SRN signal based on the flux prediction by Ando et al. |102j . 
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3. Solar flare 

Solar flares are the most energetic bursts which occur in the solar surface. Explosive release 
of energy stored in solar magnetic fields is caused by magnetic reconnections, resulting in plasma 
heating, particle accelerations, and emission of synchrotron X-rays or charged particles from the 
solar surface. In a large flare, an energy of 10^^ ergs is emitted over lO's of minutes, and the 
accelerated protons can reach energies greater than 10 GeV. Such high energy protons can produce 
pions by nuclear interactions in the solar atmosphere. Evidence of such nuclear interactions in 
the solar atmosphere are obtained via observations of solar neutrons, 2.2 MeV gamma rays from 
neutron captures on protons, nuclear de-excitation gamma rays, and possible > 100 MeV gamma 
rays from neutral pion decays. Thus, it is likely that neutrinos are also emitted by the decay of 
mesons following interactions of accelerated particles. Detection of neutrinos from a solar flare was 
first discussed in 1970's by R.Davis |104lll05j . but no significant signal has yet been found |106tll07] . 
There have been some estimates of the number of neutrinos which could be observed by large water 
Cherenkov detectors |1081 1109 J. According to [108J, about 20 neutrinos will be observed at Hyper- 
Kamiokande during a solar flare as large as the one in 20 January 2005, although the expected 
numbers have large uncertainties. Therefore, regarding solar flares our first astrophysics goal is 
to discover solar flare neutrinos with Hyper-K. This will give us important information about the 
mechanism of the particle acceleration at work in solar flares. 

4-. Indirect dark matter search 

The dark matter problem in the universe is one of the hottest topics in both particle physics and 
astrophysics. One possible solution to this conundrum is to introduce stable, weakly interacting 
massive particles (WIMPs) - which feel only the weak interaction and gravity - with sufficient mass 
to solve the structure formation problem of the universe. Therefore, many experiments are trying 
to detect WIMPs directly or indirectly all over the world. Hyper-K will be able to perform indirect 
WIMP searches as was done in Super-K |110| with improved sensitivity. In indirect searches, 
signals are from neutrinos emitted by WIMP annihilations which occur in the Sun, the Earth, or 
in the halo, where WIMPs are trapped by gravitational potentials. In this study, we focus on the 
WIMP signal from the Sun. 

These neutrinos interact with detector targets or rocks around detectors, and the secondary 
particles such as muons are detected. To distinguish those signal muons from cosmic-ray muons. 
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FIG. 65. The effective area for upward-going muon detection. Tlie red line shows Hyper-K, while the blue 
and black lines show Ice-Cube+DeepCore (at trigger level) [TTT] and KM3NET |TT2] . 

only upward-going muons are selected. The detector volume for neutrino searches using upward- 
going muons can be expressed as a neutrino effective area {A^ff), which is defined as: 

^e// = O'j^-rock X X AhK (8) 

where c7",y_rock is the cross section of neutrino-rock interaction, is the muon range, and Ahk 



is the detector area which in this case is 4550 cmx 24750 cmx2=22000 m^. Fi gure 65 shows the 
calculated effective area of Hyper-K as a function of neutrino energy. Here we reasonably assume 
the momentum threshold for muon detection in Hyper-K is similar to that of Super-K (1.7 GeV/c). 



As shown in Fig. 65, Hyper-K has a better sensitivity below a few tens of GeV than km'^-scale 
detectors thanks to its lower energy threshold for muon detection. 

In the Sun, WIMPs are scattered on hydrogen or heavy elements both with spin-dependent 
and spin-independent couplings. However, indirect searches for the annihilation signals have better 
sensitivity for the spin-dependent cross-section on hydrogen because of its dominant abundance. In 
this study, the sensitivity for WIMP signals is calculated for two cases as done in |110| 11131 1114j . In 
the first case, all WIMPs are assumed to annihilate into bb (xx ~^ ^b), the so-called "soft channel", 
while in the second case, all WIMPs are assumed to annihilate into W^W" {xx ~^ W~^W~), 
the so-called "hard channel". Since these channels give the softest and hardest neutrino spectra 
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FIG. 66. Limits on the WIMP-induced upward-going muon rate as a function of WIMP mass are shown. 
Lines show AMANDA jTTS] (green dotted Une), Ice-Cube pl4] (blue dot-dash hne), Super-K [TTO] (black 
line), and Hyper-K (red line). The shaded region is the expected flux region from DARKSUSY. In this 
analysis, all WIMPs are assumed to annihilate to bb. 



respectively, the actual sensitivity should lie somewhere between these two extremes. This approach 



is even more sensitive than the direct searches taking place on the Earth. Figures 66 and [67] show 
the experimental results of limit on WIMP-induced upward-going muon fluxes as a function of 



WIMP mass |110j and Hyper-K's expected sensitivity after 5 years of observation. In Fig. 68 



Hyper-K's sensitivity to the WIMP-proton cross section is compared to the upper limits obtained 
from both direct and indirect dark matter searches. Thank to the larger effective area of Hyper-K, 
5 years of exposure - corresponding to 90 years of exposure for Super-K - can improve the present 
sensitivity significantly. Another approach to search for WIMP annihilation neutrinos by FC events 
has been discussed elsewhere |119i 1120] . Since the neutrino energy spectrum for WIMP annihilation 
neutrinos shows a bump structure according to the WIMP mass, the reconstructed energies of FC 
events are very important to identify WIMP signals. In addition, because the atmospheric flux 
is much suppressed at higher energy (about 1/10 of at 100 GeV), the flavor information of FC 
events can be used to reduce the atmospheric neutrino background. With Hyper-K's excellent 
capabilities of energy reconstruction and neutrino flavor tagging, its soft channel WIMP-proton 
cross section sensitivity could reach 10~^^ cm^ at 10 GeV [119] . The low mass region has gotten a lot 



E Astrophysics 



97 



~ 10 



c 
o 



Hard channel 



E 

S .13 
3 

E 

Q. 
3 

■a 

(U 

u 

3 

■a 
c 



10 



10 



10 



10 



Expected flux region 
(E„ il GeV ) 



: : : AMANDA : 

— ICECUBE ; 

— SK^+ll+llj : 
rt-t-i Hy|ief-K $yr 



10 



10 



WIMP mass(GeV/c'' 



„io 



FIG. 67. Same as Fig. 66 but in this analysis, all WIMPs are assumed to annihilate to W^W 



of attention since the results of direct searches come out [115[ I12t] . Hence, Hyper-K's observation, 
which has considerably better sensitivity in this region, is likely to be very important. It should 
be noted that the larger effective area and analysis improvements in the low mass region will also 
enhance the sensitivity of searches for WIMP annihilation in the Earth or the galactic bulge as 
well as the Sun. 



5. Summary 

In this section, expectations of the Hyper-K detector for astrophysical sources are discussed. 

Quite high statistical observations of supernova neutrinos are expected for galactic supernovae, 
e.g. 170,000~260,000 events for a 10 kpc supernova. Hyper-K will give us the detailed time profile 
and temperature variation during the burst. The neutronization phase, the initial burst phase which 
emit fe's, can be detected in Hyper-K. The direction to the supernova can be determined with 
an accuracy of about 2 degrees using z/e-scattering events for a 10 kpc supernova. For supernova 
explosions outside our galaxy, we expect about 30~50 events for M31 (Andromeda Galaxy) and 
7,000~10,000 events for the Large Magellanic Cloud (LMC). If we could detect supernovae out to 
~4 Mpc, we could observe a supernova burst once every 1~3 years with high confidence. Hyper- 
Kamiokande is able to detect at least one supernova neutrino event for a supernova at 4 Mpc with 
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FIG. 68. Limits on WIMP-proton spin-dependent cross section as a function of WIMP mass. Limits from 
direct detection experiments: DAMA/LIBRA allowed region |115j (dark red and light red filled, with and 
without ion channeling, respectively), KIMS 116 (light blue dotted hne), and PICASSO |117] (light black 
dotted line) are shown. Limits from indirect detection (neutrino telescopes): AMANDA |113) (black line 
with triangles), IceCube |114j (blue line with squares), and Super-K (green line with stars, green dashed 
line) [IJjQj. Ice-Cube+DeepCore sensitivity is shown with the purple line }118j . Hyper-K sensitivity is shown 
with the red line (soft channel) and the red dashed line (hard channel). 

an efficiency of 52~69%. 

We can also investigate currently unknown properties of neutrinos by studying supernova neu- 
trinos. For example, if the neutrino mass hierarchy is inverted, the shock wave propagation may 
cause time variation of the event rate for higher energy events during the cooling phase. The energy 
dependence of arrival times at the onset of the burst would tell us the neutrino direct mass with 
an accuracy of 0.5~1.3 eV. 

Supernova relic neutrinos (SRN) will tell us the history of massive stars in the universe. In 
particular, a spectrum measurement with high statistics is necessary because the present spectrum 
of SRN is the red-shifted sum of the contribution of supernova neutrinos from every epoch of 
the universe. We expect about 310 SRN events in the energy range from 20 MeV to 30 MeV 
for 10 years of Hyper-K running, while the background from atmospheric neutrinos will be about 
2200 events without tagging neutrons. The lower energy bound (20 MeV here) is limited by 



E Astrophysics 



99 



spallation background. It is crucial to lower the detectable energy down to 10 MeV because the 
contribution from early epoch supernovae tends to be distributed at lower energy. By introducing 
0.1% gadolinium to Hyper-K, a neutron coincidence signal from the 8 MeV gamma cascade will 
remove spallation background. It will enable us to lower the energy threshold down to 10 MeV 
and also reduce atmospheric neutrino backgrounds by a factor of ~5. We would then expect about 
830 SRN events per 10 years in E=10-30MeV assuming 67% efficiency for tagging neutrons. 

For the solar flare neutrinos, even though there are large uncertainties in predictions, the existing 
observations of accelerated particles and secondary particles support the supposition that neutrinos 
are also produced during a solar flare. The first detection of solar flare neutrinos by Hyper-K will 
be a strong test of neutrino emission models. For the indirect dark matter search, 5 years of 
observation by Hyper-K will yield a sensitivity to the WIMP-proton spin dependent cross section 
below 10~^^ cm^ around 10 GeV for the soft channel, and 10~^° cm^ around 100 GeV for the 
hard channel. It can improve the present sensitivity significantly for lighter WIMPs below ^■lOO 
GeV. Other astronomical neutrino searches such as GRB neutrinos, galactic diffuse neutrinos, etc. 
can be also performed in Hyper-K with its large acceptance in a wide energy range as well as its 
excellent performance of event reconstruction. 
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F. Neutrino geophysics 

In this section, the radiographic measurements of the Earth with Hyper-Kamiokande will be 
discussed. 

1. Geophysical motivation 

The Earth's internal structure and chemical composition have been estimated by analyzing 
seismic data to derive seismic velocities inside the Earth, in conjunction with many laboratory 
experiments and model calculations. Until recent years, these seismic waves had been the only 
probe that could penetrate the Earth. In order to construct the standard Earth model a number 
of models have been considered. PREM [69] and 1066A |122j are two major models, both of 
which describe the spherically symmetric structure of the Earth. PREM and 1066A are both 
parametrized models where the seismic velocities Vp{x) and Vs{x) are assumed to have a certain 
relationship with the density p{x) via the gravity g, bulk modulus k, and so on, and thus there 
is an uncertainty regarding the absolute density distribution. It is therefore critically important 
to find an independent method to directly measure the density inside the Earth. The question is: 
what kind of method can be used for this purpose? Drilling and core sampling enable us to directly 
examine material inside the Earth. However, considering the fact that the world's deepest such 
sample is 12 km deep, it is likely impossible, not to mention expensive and dangerous, to attempt 
to reach the core (deeper than 3000 km) by this method. 

Recent observation of anti-neutrinos from the decay of radiogenic isotopes inside the Earth by 
the KamLAND jl23t I124j and Borexino |125j experiments opened novel possibilities to investi- 
gate Earth's interior using neutrinos. In what follows, another approach of neutrino geophysics 
using atmospheric neutrinos detected by Hyper-K is discussed. The observation of neutrinos and 
measurement of their absorption or differential oscillations in the Earth will provide unique infor- 
mation about its density structure. Furthermore, the result will be more easily interpreted than 
the conventional seismological method which has intrinsic uncertainties. 

We estimate the accuracy necessary to provide useful information on geophysics, taking as an 
example measuring the density of the Earth's core. The main component of the Earth's core is 
assumed to be metal iron. Seismic measurements indicate that the density of the core is 10% less 
than that of pure iron |126) . On the other hand, siderophile elements, elements that can easily 
dissolve into iron, are found to have extremely low abundance in the mantle. Ringwood |127| 
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explained the reason why the amount of the siderophile elements is extremely low by using his 
model of element transportation between the mantle and the core. The siderophile elements that 
are lighter than iron can dissolve into the core iron, and as a result their abundance in the mantle 
decreases. This seemed to give an explanation of the core density being lighter than pure iron. 
However, when Ringwood surveyed the amount of the siderophile elements dissolved into the core 
it was found that the core abundance of the siderophile elements could not explain its anomalous 
density. The low density core remains, despite active inquiries and lively discussions, a mystery. 

Thus, a measurement of the Earth's density with the accuracy of 10% or better can provide 
useful information on geophysics. If the accuracy is further improved, other topics in geophysics 
can be pursued. Following are examples of topics of great interest: 

1. the density difference between the inner core and the outer core: this difference reflects the 
partition coefficient of light elements between the inner and outer cores, which is important 
to understand the energy source of the geodynamo. A measurement accuracy of 5% will be 
required for this subject. 

2. the value of the Cq parameter of Roche's law: this parameter contains the history of the 
inner core [128j . The subject of interest is the convection of the inner core. If this value 
is precisely determined, the gradient of chemical composition can be discussed. If there is 
convection, the chemical composition would be homogeneous within the core, while if not it 
would show a gradient distribution. 

3. whether the inner core is spherically or cylindrically symmetric: the seismic structure indi- 
cates that the inner core is cylindrically symmetric rather than spherically. Also, there is 
a report that indicates asymmetry between the east and the west hemispheres of the inner 
core, although it is difficult to quantitatively evaluate this difference. 

2. Current studies 

Neutrino radiography for surveying the internal structure of the Earth was first discussed more 
than 25 years ago [12911131] . and many ideas have been proposed since then. However, because 
those ideas are based on gigantic accelerators or use of (thus far unobserved) galactic neutrinos, 
they are difficult to realize. Recently, there has been an attempt to use atmospheric neutrinos for 
neutrino radiography, based on a simulation of atmospheric neutrino events that can be collected 
with the IceCube neutrino detector. It is found that the density difference between the core and the 
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mantle can be determined with an accuracy of 20% in 10 years |132) . Following this calculation, the 
data taken with IceCube was analyzed. Based on the one year data taken with 40 strings (about 
half of the full IceCube detector), the density difference between the core and the mantle can be 
confirmed at the la CL |133j . Based on a simulation of neutrino radiography with KMSNeT |134] . 
it was reported that the average density of the core can be determined with accuracy of 6% {la 
CL) in 10 years. 



3. Radiography with Hyper- Kamiokande 

The absorption neutrino radiography utilizes very high energy neutrinos with energies above 
10 TeV, and thus a very large detector is naturally required. The IceCube and the KMSNeT 
detectors are designed specifically to detect very high energy neutrinos. They are not optimized 
for detection of low energy neutrinos because the optical modules (or photo sensors) are placed too 
sparsely. 

Although Hyper-Kamiokande does not have such a large volume as IceCube or the KMSNeT, its 
capability to detect a large number of (lower energy) atmospheric neutrinos enables us to perform 



the neutrino oscillation radiography using the MSW effect [65^16 7j . Figure 69 shows the ratio of 
atmospheric Ve flux with MSW effect to that without MSW effect as a function of zenith angle and 
neutrino energy. The largest change of the atmospheric neutrino oscillation probability by the MSW 
effect in the Earth's core can be seen in the energy range between 5 and 10 GeV as discussed in 



Sec. Ill B The density resolution is found to depend strongly on the neutrino oscillation parameters 
^23 J Am32, ^13, and 5. The flux of neutrinos in the energy range between 5-10 GeV is more than 
10^ times greater than the flux of such very high energy neutrinos (assuming the differential cosmic 
ray spectral index of 2.7). Given the precisely known oscillation parameters, Hyper-K would be 
able to perform the neutrino oscillation radiography with considerably higher statistics than these 
km^-scale detectors based on MSW effect. 
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FIG. 69. Ratio of atmospheric i^e flux with MSW effect to that without MSW effect, as a function of neutrino 
energy and zenith angle, sin^ ^?23 = 0.5, sin^ 2^13 = 0.1, and normal mass hierarchy are assumed. Other 
parameters are listed in Table |XI[ 
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